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ABSTRACT 

We present a new model to describe the galaxy-dark matter connection across cosmic time, which 
unlike the popular subhalo abundance matching technique is self-consistent in that it takes account of 
the facts that (i) subhalos are accreted at different times, and (ii) the properties of satellite galaxies 
may evolve after accretion. Using observations of galaxy stellar mass functions out to z ~ 4, the 
conditional stellar mass function at z ^ 0.1 obtained from SDSS galaxy group catalogues, and the 
two-point correlation function (2PCF) of galaxies at z ~ 0.1 as function of stellar mass, we constrain 
the relation between galaxies and dark matter halos over the entire cosmic history from z ~ 4 to the 
present. This relation is then used to predict the median assembly histories of different stellar mass 
components within dark matter halos (central galaxies, satellite galaxies, and halo stars). We also 
make predictions for the 2PCFs of high-z galaxies as function of stellar mass. Our main findings are 
the following: (i) Our model reasonably fits all data within the observational uncertainties, indicating 
that the ACDM concordance cosmology is consistent with a wide variety of data regarding the galaxy 
population across cosmic time, (ii) At low-z, the stellar mass of central galaxies increases with halo 
mass as M°-^ and at the massive and low-mass ends, respectively. The ratio M^_c/M reveals a 

maximum of ^ 0.03 at a halo mass M ^ lO^^'* /i~^M0, much lower than the universal baryon fraction 
(~ 0.17). At higher redshifts the maximum in M^^^M remains close to ~ 0.03, but shifts to higher 
halo mass, (iii) The inferred time-scale for the disruption of satellite galaxies is about the same as 
the dynamical friction time scale of their subhalos. (iv) The stellar mass assembly history of central 
galaxies is completely decoupled from the assembly history of its host halo; the ratio M^,^c/M initially 
increases rapidly with time until the halo mass reaches - W^h~^MQ, at which point M^^c/M ^ 0.03. 
Once M > 10^^ /i~^Mq, there is little growth in M*_c, causing the ratio M^^c/M to decline. In Milky- 
Way sized halos more than half of the central stellar mass is assembled at z < 0.5. (v) In low mass 
halos, the accretion of satellite galaxies contributes little to the formation of their central galaxies, 
indicating that most of their stars must have formed in situ. In massive halos more than half of the 
stellar mass of the central galaxy has to be formed in situ, and the accretion of satellites can only 
become significant at z < 2. (vi) The total mass in halo stars is more than twice that of the central 
galaxy in massive halos, but less than 10 percent of M*_c in Milky Way sized halos. (vii) The 2PCFs 
of galaxies on small scales holds important information regarding the evolution of satellite galaxies, 
and at high-z is predicted to be much steeper than at low-z, especially for more massive galaxies. 
We discuss various implications of our findings regarding the formation and evolution of galaxies in a 
ACDM cosmology. 

Subject headings: dark matter - large-scale structure of universe - galaxies: halos 



1. INTRODUCTION 

In recent years, much effort has been made to estab- 
lish the statistical connection between galaxies and dark 
matter halos, as parameterized via the conditional lu- 
minosity function (CLE) (Yang, Mo & van den Bosch 
2003) or the halo occupation distribution (HOD) (Jing 
et al. 1998; Peacock & Smith 2000). This galaxy-dark 
matter connection describes how galaxies with different 
properties occupy halos of different mass, and yields im- 
portant insight into how galaxies form and evolve in 
dark matter halos. In practice, the various methods 
that have been used to constrain the galaxy-dark mat- 
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ter connection (galaxy clustering, galaxy-galaxy lensing, 
galaxy group catalogues, abundance matching, satellite 
kinematics) use the fact that the halo properties, such 
as mass function, mass profile, and clustering, are well 
understood in the current ACDM model of structure for- 
mation (e.g.. Mo, van den Bosch & White 2010). 

At low-redshift, large redshift surveys, such as the two- 
degree Field Galaxy Redshift Survey (2dFGRS; CoUess 
et al. 2001) and the Sloan Digital Sky Survey (SDSS; 
York et al. 2000), have provided accurate estimates of 
the luminosity and stellar mass functions of galaxies (e.g., 
Norberg et al. 2002b; Blanton et al. 2003; Li & White 
2009), of their two-point correlation functions (2PCF) 
as function of various galaxy properties (e.g., Norberg 
et al. 2002a; Zehavi et al. 2005, 2011; Wang et al. 2007), 
their satellite kinematics (e.g., van den Bosch et al. 2004; 
More et al. 2009, 2011), and even of their excess sur- 
face densities, a measure for the tangential shear caused 
by gravitational lensing due to their mass distributions 
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(e.g., Mandelbaum et al. 2005). All these results have 
been used to infer how galaxies with different properties 
are distributed in halos of different masses (e.g. Jing et 
al. 1998; Peacock & Smith 2000; Yang et al. 2003; van 
den Bosch et al. 2003a, 2007; Zheng et al. 2005; Tinker 
et al. 2005; Mandelbaum et al. 2006; Brown et al. 2008; 
More et al. 2009, 2011; Caceiato et al. 2009; Neistein et 
al. 2011a, b; Avila-Reese & Firmani 2011). In addition, 
these large galaxy redshift surveys can also be used to 
identify galaxy systems (groups), defined as those galax- 
ies that share a common dark matter host halo (Yang 
et al. 2005a; 2007). Such group catalogs can be used 
to examine the halo - galaxy connection even more di- 
rectly than the methods mentioned above (e.g.. Eke et 
al. 2004; Yang et al. 2005b, 2008, 2009b). These analyses 
have revealed a number of important properties regard- 
ing the relation between galaxies and their dark matter 
halos. The stellar mass to halo mass ratio has a mini- 
mum for halos of ~ 10"'^^ /i~^M0, and increases rapidly 
towards both lower and higher halo masses (Yang et al. 
2003, 2005b, 2008, 2009b; van den Bosch et al. 2003a; 
Tinker et al. 2005; Leauthaud et al. 2012), suggesting 
that galaxy formation is most efficient in halos with a 
present-day mass of ^ lO'^^ h~^MQ, i.e., these halos have 
the highest integrated star formation efficiencies. Yet, 
their total stellar masses are only a few percent of that 
of the halo, indicating that the overall star formation 
efficiency is very low. 

Ideally, we would like to carry out similar analyses at 
various redshifts, so as to investigate how star forma- 
tion proceeds and how galaxies assemble as their host 
halos grow within the cosmic density field. Such analy- 
ses are now becoming possible at intermediate redshift, 
2; ~ 1, where reliable luminosity and stellar mass func- 
tions of galaxies have been obtained from various redshift 
surveys, such as the DEEP2 survey (Davis et al. 2003), 
the COMBO-17 survey (Wolf et al. 2004), VVDS (Le 
Fevre et al. 2005), and zCOSMOS (Lilly et al. 2007). 
In addition, these surveys have been used to measure 
the 2PCFs of galaxies as function of their luminosity, 
stellar mass and/or color (e.g., Daddi et al. 2003; Coil 
et al. 2006; Phleps et al. 2006; Polio et al. 2006; Mc- 
Cracken et al. 2008; Meneux et al. 2008, 2009; Foucaud 
et al. 2010; de la Torre et al. 2010). These observa- 
tions have prompted a series of investigations into the 
galaxy-dark matter connection and its evolution between 
2; ~ 1 and the present (e.g.. Bullock et al. 2002; Mous- 
takas & Somerville 2002; Yan et al. 2003; Zheng 2004; 
Lee et al. 2006; Hamana et al. 2006; Cooray 2005, 2006; 
Cooray & Ouchi 2006; Conroy et al. 2005, 2007; White 
et al. 2007; Zheng et al. 2007; Conroy k Wechsler 2009; 
Wang & Jing 2010; Wetzel & White 2010; Wang & Jing 
2010; Leauthaud et al. 2011; Wake et al. 2011). 

With the advent of deep, miilti wave-band surveys, 
it has even become possible to estimate the luminos- 
ity/stellar mass functions of galaxies out to z ^ 8 (e.g. 
Drory et al. 2005; Fontana et al. 2006; Perez-Gonzalez 
et al. 2008; Marchesini et al. 2009; Stark et al. 2009; 
Bouwens et al. 2011). However, the data samples are 
still small (and hence subject to cosmic variance), with 
large discrepancies among different measurements (see 
Marchesini et al. 2009). Furthermore, since reliable clus- 
tering measurements are in general unavailable for these 
high redshift galaxy samples, it is not possible to carry 



out the same HOD/CLF analyses for these high-z galax- 
ies as for galaxies at low z. Nevertheless, attempts have 
been made to establish the relation between galaxies and 
their dark matter halos out to high z using a technique 
known as abundance matching, in which galaxies of a 
given luminosity or stellar mass are linked to dark mat- 
ter halos of given mass by matching the observed abun- 
dance of the galaxies in question to the halo abundance 
obtained from the halo mass function (typically also ac- 
counting for subhalos). This approach was first used by 
Mo & Fukugita (1996) and Mo, Mao & White (1999) 
to model the number density and clustering of Lyman- 
break galaxies. More recently, several studies have used 
abundance matching techniques to probe the galaxy-dark 
matter connection out to z ~ 5 (e.g.. Vale & Ostriker 
2004, 2006; Conroy, Wechsler & Kravtsov 2006; Shankar 
et al. 2006; Conroy & Wechsler 2009; Moster et al. 2010; 
Guo et al. 2010; Behroozi, Conroy & Weschler 2010). 

An important aspect of abundance matching is the 
treatment of dark matter subhalos. It is usually as- 
sumed that a central galaxy resides at the center of each 
halo, and orbiting around it are satellite galaxies associ- 
ated with the subhalos of the host halo. These subhalos 
were distinctive (host) halos themselves before they were 
accreted into their hosts, namely satellite galaxies were 
themselves central galaxies before their host halo became 
a subhalo. The modeling of the total galaxy population 
can therefore be separated into two parts: (i) the forma- 
tion of central galaxies in dark matter halos at different 
redshifts, and (ii) the accretion and evolution of satel- 
lite galaxies in their host halos. With the use of high- 
resolution numerical simulations and extended Press- 
Schechter theory (Bond et al. 1991; Lacey & Cole 1993), 
the properties of the subhalo population have now been 
determined with great accuracy (e.g. Gao et al. 2004; De 
Lucia et al. 2004; Tormen et al. 2004; van den Bosch et 
al. 2005a; Weller et al. 2005; Diemand et al. 2007; Gio- 
coU et al. 2008; Li & Mo 2009). If satellite galaxies are 
indeed associated with subhalos, these properties should 
be relevant to the modeling of the galaxy population. In 
particular, since subhalos were host halos before accre- 
tion, the subhalo property that seems most relevant is 
its mass at the time of accretion. Consequently, many 
abundance matching studies to date have linked galaxies 
of property V (i.e., galaxy luminosity or stellar mass) to 
dark matter (sub)-halos of mass M using 

roo poo 

I n^{V,z)dV= / n^{M,z)dM 
Jv Jm 

/.oo 

+ / risub (ma, 2) drria. (1) 
Jm 

Here ng{'P, z) is the comoving number density of galaxies 
of property V at redshift z, nh{M,z) is the halo mass 
function at redshift z, and Usuhi'nT'a, z) is the comoving 
number density of subhaloes identified at redshift z which 
at accretion have a mass ma- Hence, in such abundance 
matching the property of a satellite galaxy is assumed to 
depend only on its halo mass at accretion. 

Although subhalo abundance matching yields galaxy 
correlation functions that are in remarkably good agree- 
ment with observations (e.g., Conroy et al. 2006; Guo 
et al. 2010; Wang & Jing 2010), it implies a particu- 
lar path for the evolution of satellite galaxies. Indeed, 
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assuming that the stellar masses of satellite galaxies de- 
pend only on their halo mass at accretion implies either 
that the relation between V and halo mass M is inde- 
pendent of the time when the subhalo is accreted, and 
the evolution after accretion is independent of the host 
halo into which the subhalo has been accreted and how 
long ago it has been accreted, or that the effects of differ- 
ent accretion times and subsequent evolutions in different 
hosts conspire to give a stellar mass that depends only 
on the mass of the subhalo at accretion. Such evolution 
of satellites is neither physically motivated nor the only 
possible path which satellite galaxies could take. In fact, 
applying the abundance matching technique as described 
above to data at different z results in a P{M) relation 
that changes with redshift (e.g., Conroy et al. 2006; Con- 
roy & Wechsler 2009; Moster et al. 2010; Wang & Jing 
2010). If one accepts that V for central galaxies depends 
on both M and z, then the properties of satellite galaxies 
do not just depend on the host halo mass at accretion, 
TTia, but also on the accretion redshift, Za- Furthermore, 
after accretion a satellite galaxy may lose or gain stellar 
mass and even be disrupted due to tidal stripping and 
disruption, such that V of satellite galaxies may depend 
on other properties in addition to both to^ and Za- So 
far these effects have not been modeled in any detail in 
abundance-matching studies. 

In this paper, we develop a new, self-consistent ap- 
proach to model the relation between galaxies and dark 
matter halos over cosmic time. Our approach is based 
on the model of Yang et al. (2011), which gives subhalo 
abundance as a function of both mass at accretion and 
accretion time. This allows us to use a galaxy-halo rela- 
tion that depends on subhalo accretion time and to follow 
the dynamical evolution after their accretion. Further- 
more, since satellites observed in halos (galaxy systems) 
at the present time are 'fossil records' of both the for- 
mation of central galaxies at high-^ and the subsequent 
dynamical evolution, the observed abundance and clus- 
tering of satellite galaxies at the present time can be used 
to constrain the galaxy - halo relations at different red- 
shifts. In the present work, we use our model together 
with various observational data to obtain the stellar mass 
- halo mass relation as a function of redshift from z ^ 4 
to the present. This paper is organized as follows. In Sec- 
tion [2] wc outline the halo model ingredients to be used. 
In Section [3] we describe how we model the stellar mass 
function, the conditional stellar mass function, and the 
two point correlation function of galaxies. The observa- 
tional constraints used for our analysis are presented in 
Section S) The results are presented in Section [SJ while 
Section |6] describes the implications regarding the assem- 
bly of galaxies and the evolution of the galaxy dark mat- 
ter connection. We summarize our results in Section [71 
which also includes a detailed discussion regarding vari- 
ous implications for galaxy formation and evolution. 

Throughout this paper we will mostly focus on 
a ACDM 'concordance' cosmology whose parameters 
are consistent with the seventh-year data release of 
the WMAP mission (Komatsu et al. 2011; hereafter 
WMAP7). However, in order to investigate the pos- 
sible cosmology dependence of our results, we also 
present results for the following set of ACDM cosmolo- 
gies: WMAPl, with parameters consistent with the first- 
year data release of the WMAP mission (Spergel et 



TABLE 1 

ACDM COSMOLOGICAL MODELS USED IN THIS PAPER. 



Name 






71s 


h 




WMAPl 


0.30 


0.70 


1.00 


0.70 


0.90 


WMAP3 


0.238 


0.762 


0.951 


0.73 


0.75 


WMAP5 


0.258 


0.742 


0.963 


0.719 


0.796 


WMAP7 


0.275 


0.725 


0.968 


0.702 


0.816 


Millennium 


0.25 


0.75 


1.00 


0.73 


0.90 



al. 2003); WMAP3, with parameters given in Spergel et 
al. (2007); WMAP5, with parameters given in Dunkley et 
al. (2009), and finally a cosmology ('Millenium') with pa- 
rameters that are identical to those adopted in Springel 
et al. (2005) for the Millennium Simulation. Table [T] lists 
the parameters for all these cosmologies. 

2. DARK MATTER HALOS AND SUBHALOS 

The main goal of this paper is to obtain a self- 
consistent model for the link between galaxies and dark 
matter halos accross cosmic time. In this section we de- 
scribe our model ingredients for halos and subhalos that 
are needed for our investigation. 

2.1. Dark Matter Halo Mass Function 

The mass function of dark matter halos, nh(M, z) dM, 
describes the number density of dark matter halos of 
mass M at redshift z. The Press-Schechter formalism 
(Press & Schechter 1974) yields an analytical estimate for 
nh(M, z), and we use the revised form given in Sheth, Mo 
& Tormen (2001) based on the ellipsoidal collapse model: 



n^{M,z)dM^^vf{u) 



diner 



dlnM 



dM, 



(2) 



where p is the mean matter density of the Universe at 
z, cr(M) is the mass variance, and v = 6c{z) / a{M) with 
5c{z) the critical over-density required for collapse at red- 
shift z. The function f{v) is given by 

-/M-2Afl + f^V'^cxpf-^) (3) 



with v' ^ y/^v, a = 0.707, q = 0.3 and A « 0.322. The 
resulting mass function has been shown to be in excellent 
agreement with numerical simulations, as long as halo 
masses are defined as the masses inside a sphere with an 
average over-density of about 180 that of the background 
mass density (Sheth & Tormen 1999; Jenkins et al. 2001). 
As an illustration, the left-hand panel of Fig. [T] shows the 
halo mass functions at several redshifts for the WMAP7 
cosmology, predicted using Eq. ([2|). 

2.2. Halo Density Profile 

Throughout this paper we assume that dark matter 
halos are spherical and follow an NEW density profile 
(e.g., Navarro, Frenk & White 1997): 



p{r) 



5-p 



(r/r,)(l + r/r,)2' 



(4) 



where is a characteristic radius, and (5 is a dimension- 
less amplitude which can be expressed in terms of the 
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Fig. 1. — Left panel: the halo mass functions predicted using the Sheth Mo & Tormen (2001) model for the WMAP7 cosmology at 
different redshifts, as indicated. Right panel: the bias parameter of dark matter halos, again predicted using the Sheth, Mo & Tormen 
(2001) model for the WMAP7 cosmology at the same redshifts as those indicated in the left panel. 



halo concentration parameter c = rigo/rg as 



(5 = 



180 



3 ln(l + c) - c/(l + c) 



(5) 



Numerical simulations have shown that c is correlated 
with halo mass (e.g., NFW; Eke et al. 2001; Jing & 
Suto 2000; Bullock et al. 2001; Zhao et al. 2003; 2009; 
Maccio et al. 2007). Throughout this paper we use the 
concentration- mass relation of Zhao et al. (2009), prop- 
erly corrected for our definition of halo mass. 

2.3. Halo Bias 

Dark matter halos are biased tracers of the dark matter 
mass distribution. The amplitude of this bias depends 
on halo mass and is expressed via the halo bias function, 
hi,{M) (e.g.. Mo & White 1996, 2002). This allows one 
to write the two point correlation function between halos 
of masses Afi and M2 on large scales as 



ehh(r|Afi, A/2) = 6h(Mi) 6i,(A/2) edm(r) 



(6) 



where Cdm(?') is the two-point correlation function of the 
(non-linear) dark matter mass distribution. Throughout 
this paper we adopt the halo bias function of (Sheth, Mo 
& Tormen 2001), which is given by 



bias factor, bh{M) at several redshifts for the WMAP7 
cosmology, predicted using Eq. ((T]). 

2.4. The Suhhalo Population 

An important ingredient of our model describing the 
link between galaxies and dark matter halos is a statis- 
tical description for the population of dark matter sub- 
halos. In particular, we need to know the distribution 
of accretion masses, ma, and accretion redshifts, Za, for 
the population of subhalos, as a function of the mass 
of its host halo. Let nsub(™a, Za\M, z) dm^ dln(l -I- Za) 
denote the number of subhalos in a host halo of mass 
M at redshift z, as a function of their accretion masses, 
nia, and accretion redshifts, Zq. We will refer to ngub as 
the subhalo mass function (SHMF). In the past, the sub- 
halo mass function has typically been studied using A^- 
body simulations or Monte-Carlo realizations of the ex- 
tended Press-Schechter formalism (e.g., Sheth & Lemson 
1999; Somerville & Kolatt 1999; Cole et al. 2000; van den 
Bosch et al. 2005a; Giocoli et al. 2008a; Cole et al. 2008; 
Parkinson et al. 2008; Fakhouri & Ma 2008; Fakhouri 
et al. 2010). Recently, however, Yang et al. (2011; 
hereafter Yll) presented a fully analytical model for 
nsnhifna, Za\M, z) which is the one we will adopt here. 
Following their notation we have that 



&h(M) = l 



^(ai/2) 
(ai/2)= 



nsnhfi{^a,Za\M,z) = Na{s a, 5a\S , 5) 



(9) 



(ai/2)'=-h6(l-c)(l-c/2) 



(7) 



with a = 0.707, h = 0.5, c = 0.6 and v = 8^1 a{M). Note 
that our halo mass function ([3]) and halo bias function 
([7|) obey the normalization condition 



(8) 



which expresses that the distribution of dark matter is, 
by definition, unbiased with respect to itself. As an il- 
lustration, the right-hand panel of Fig. [1] shows the halo 



where Sa = a'^ima), S = (T^(Af), Sa = Sc{za), 5 = S^z), 
and the subscript '0' on the SHMF indicates that this 
function ignores higher order subhalos (subhalos of sub- 
halos, etc; see below). Yll considered three different 
models for J\fa{sa,Sa\S,d), and we adopt their Model 
III which is the most accurate. For host halos of given 
(Af, z), this model uses the mean halo assembly history 
given by the fitting formula of Zhao et al. (2009) and a 
log-normal model to describe the scatter among different 
halos (see Section 3 of Yll for details). Note that this 
SHMF only includes subhalos that are directly accreted 
onto the main branch of the host halo (hence the sub- 
script '0'). However, since the SHMF given by Eq. ^ 
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also applies to subhalos before their accretion (when they 
themselves were host halos), we can in principle calcu- 
late the SHMF including higher level sub-halos, i.e., sub- 
halos within sub-halos, etc. For example, the SHMF of 
sub'-subhalos (ith level subhalos) can be written as 



nsuh,i{ma,Za\M,z) = J J nsnh,i-lima,Za\mi^i,Zi.i) 

xnsub,o(™i-i, Zj-i|M, z) dmj_i dln(l + . (10) 

Thus the total subhalo population can be described using 



nsuh{ma,Za\M,z) 



nsnh.i{'>^a,Za\M,z) 

i=0 



(11) 



In practice, we find that the sub*-subhalos for i > 2 
contribute less that 5% to ?T.sub(wa, Za\M, z) in the mass 
range of interest here (> 10^° /i~^Mq), and we therefore 
adopt A^max = 1 in what follows (see also Yang, Mo & 
van den Bosch 2009a). 

3. MODELING MASS FUNCTIONS AND 
CORRELATION FUNCTIONS 

3.1. The Stellar Mass Functions 

In the conditional luminosity (stellar mass) function 
model developed by Yang et al. (2003), the stellar mass 
function of galaxies at a given redshift z can be written 

as 



$(Af„z) 



$(M,|M,z)rih(M,z)dM, (12) 



where $(M*|M, z) is the conditional stellar mass func- 
tion (CSMF), which gives the mean number of galaxies 
with stellar masses in the range ± dMH</2 hosted by 
halos of mass M at z. In general, we split the CSMF in 
two parts, 

$(M,|M, z) = $c(M*|Af, z) + $s(to*|A^, z) . (13) 

where i'c(-Af*|M, z) and $s(™*|-Af, z) are the contribu- 
tions from the central and satellite galaxies, respectively. 
Once $(M»|M, z) is known, the average halo occupation 
number (HON) of galaxies within a given stellar mass 
range, M*_i < < M,^2 can be written as 



{N\M,z) 



M, : 



$(M,|M,z)dM, 



(14) 



Using Eq. (|13p we write 

{N\M, z) = {N,\M, z) + {N,\M, z) , (15) 

where Nc and Ns are the occupation numbers of central 
and satellite galaxies, respectively. 

3.1.1. The Conditional Stellar Mass Function of Central 

Calaxies 

Based on the results of Yang et al. (2009b, hereafter 
Y09b), we assume that the CSMF of central galaxies is 
given by a lognormal distribution: 



$e(M*|M,z) = 



1 



27r(T, 



■exp 



(logAf,/M,,e)' 



2a2 



(16) 



where log Af,.c is the expectation value of the (10-based) 
logarithm of the stellar mass of the central galaxy and CTc 



is the dispersion. For simplicity, throughout the paper, 
we assume Uc to be independent of halo mass. 

Numerous studies have attempted to constrain the re- 
lation between Af,.c and M and its evolution with red- 
shift (see Section [1}. Following Y09b, we assume that 
this relation has the form of a broken power-law: 



(Af/Af] 



\a+l3 



{1 + M/Mi)f> 



(17) 



so that AU^c oc M°'+^ (A/,,c oc Af") for M < Mi 
(Af ^ Ml). This model contains four free parameters: 
an amplitude M^^, a characteristic halo mass. Mi, and 
two power-law slopes, a and (3. Note that all four param- 
eters may depend on redshift, which we will parameterize 
accordingly in Section below. 

3.1.2. The Conditional Stellar Mass Function of Satellite 

Galaxies 

The conditional stellar mass function for satellite 
galaxies can formally be written as, 

$s(w*|Ar,z) 

oo oo oo Ml 

= y d log m'^ J dma J ^ J dMa J dr] 

2 

P{m^\m'^,ma; Za, Ma] ri, z) 

>I>e(m',|ma, Za, z) n,sub("la, Za\M,z) 
P{Ma,Za\M,z)P{r]). (18) 

Here risuhi'ma, Za\M, z) is the SHMF described in Sec- 
tion 12 .4( ^oim'Jma, Za, z) is the CSMF of galaxies 
that were accreted at redshift Zq, taking into account 
their subsequent mass evolution to the redshift z in 
question due to star formation and stellar evolution; 
P{m^\m'^, Za', Ma', rria', rj, z) is the probability that a satel- 
lite galaxy whose stellar mass would be m'^ in the absence 
of stripping or disruption ends up with a stellar mass to* 
at redshift z due to such environmwental effects. This 
last probability is written in a form that depends explic- 
itly on the masses of the subhalo and host halo at the 
accretion epoch (toq and Ma, respectively) as well as on 
the orbital circularity of the subhalo, -q (to be defined 
below). Note that the integrand includes the probabil- 
ity distributions P{'q) and P(A/q, Za| Af, z). The latter 
describes the probability that the main progenitor of a 
halo of mass M at redshift z has a mass Ma at red- 
shift Za > z; hence P{Ma, Za\M, z) describes the mass 
assembly history of the host halo. Note also that the 
integration over Za implies that the CSMF of satellite 
galaxies at redshift z depends on the CSMF of centrals 
at all redshifts Za > z. The orbital circularity is defined 
as rj = i / ic{E) {0 < rj < 1). Here j is the orbital spe- 
cific angular momentum and jc{E) is the specific angular 
momentum of a circular orbit that has the same orbital 
energy, E, as the orbit in question. Numerical simula- 
tions have show that the probability distribution of rj of 
dark matter subhaloes is well approximated by 



P(77)cxr;i-2(l-,y)i 



(19) 



independent of redshift or halo mass (e.g., Zentner et 
al. 2005). 
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An important quantity that must somehow enter in a 
description of -P(m*|m'^, iria; Za', Ma; i], z) is the dynami- 
cal friction time scale, ^df , defined as the time interval be- 
tween the accretion of a subhalo and the epoch at which 
it is either tidally disrupted or cannibalized by the cen- 
tral galaxy. Using A^-body simulations, Boylan-Kolchin 
et al. (2008) have shown that 



Here 



= 0.216- ^ ' ' 



^dyn 



ln(l + Ma/nia) 



,1.9r, 



(20) 



which is the functional form we adopt throughout. Here 
Tviriza) is the virial radius of the host halo at Za, and 
Tdyn ~ 0.1-ff~^(2;a) is the dynamical time of the halo at 
Za- As an approximation, we assume 



rdE) = r^ir{Za) , 



(21) 



so that td{ = tdf{ma/Ma, Za,ri), and can thus be evalu- 
ated for the integrand in Eq. (fTS)) . 

In this paper we assume that a galaxy after becoming 
a satellite can gain stellar mass due to star formation 
and suffer mass loss due to passive evolution. As a sim- 
ple model, we assume the overall evolutionary effect can 
be modeled as a function of the masses m*^z and m*_a, 
where to*_z is the expected median stellar mass of central 
galaxies in halos of mass ma at redshift z. Specifically, 
we write the evolution of the median stellar mass of a 
satellite galaxy as 



TO. 



(1 



cm. 



(22) 



where c is a parameter which may depend on z. Thus, 
if c = 0, then m'^ = m^^a so that the stellar mass of a 
satellite is equal to its original mass at accretion. On the 
other hand if c = 1, then m'^ = m^^z so that the stellar 
mass of a satellite is the same as that of central galaxy 
of the same halo mass at the redshift, z, in question. 
The c = 1 case corresponds to the assumption adopted 
in the conventional abundance matching represented by 
equation ([1]). By treating c as a free parameter, we can 
examine how the evolution may deviate from the simpli- 
fied assumptions. Given this assumption, we write 



^e{m'Jma,Za,z) = 



1 



2n(T' 



■cxp 



(logm'Jm'J^ 



2a'2 



(23) 

As a simple model, we assume a'^ to be the same as CTc 
that specifies the CSMF of central galaxies (see Section 
3.3). Thus, in the special case of c = 0, ^e{m'Jma, Za, z) 
has the same form as the CSMF of central galaxies at 
the accretion redshift z^. In general, $0 is a lognormal 
distribution with a median rn'^ and a dispersion given by 

For the stripping and disruption of the satel- 
lite galaxies, we adopt a simple form for 
P{m^,\m^,^cm, Za;ma; Ma;rj, z), assuming that signif- 
icant tidal stripping occurs only for the dark matter 
subhalo, but not for the stellar component that repre- 
sents the satellite galaxy, at least not until it is either 
destroyed by the tidal forces of, or cannibalized by, the 
central galaxy. In this case, we have that 



P{m^\m'^,Za;ma;Ma;r],z) = 

(5^(m, — m',) 




if At < pt tdf /24) 
otherwise . ^ ' 



At = t{z) - t{Za) ., 



(25) 



is the time interval between z and Za (i.e., the time since 
accretion). Thus we have two free parameters c and pt 
that we seek to constrain using the observed satellite 
population. The parameter pt is intended to account for 
the tidal stripping effect. The extreme case of pt = oo 
corresponds to no tidal stripping/disruption of satellite 
galaxies, while the case of pt = corresponds to instan- 
taneous disruption of all the satellite galaxies. Under the 
above assumptions, Eq. (fT8| reduces to 



M 



$s(m*|M, z) 



dnir, 



dZa 



M 



1 



dMa / d77 



z 

$o(TO*|TOa, Za, z) r7,sub(TOa, Za\M, z) 

P{Ma, Za\M, z) P{r^) Q{pt tdf - At) 



(26) 



with Q[x) the Heaviside step function, which is straight- 
forward to compute numerically. 

For the zth-order sub-subhalos (note that we only con- 
sider i = 1), we assume that the dynamical friction time 
is also given by equation (PO)) . but with Ma referring to 
the mass of the main branch of the (m^^i, Zi^i) halo at 
the accretion time Za = Zi. This dynamical friction time 
has to be compared with At = t{zi^i) — t{zi) to decide 
if the {rrii, Zi) halo and its central galaxy were disrupted 
already in the (i — l)th-order subhalo. Their subsequent 
disruption then follows that of the [i — l)th-order (here 
Oth-order) subhalos. 

3.2. The Two-Point Correlation Function of Galaxies 

In modeling the two-point correlation function (2PCF) 
of galaxies we follow the procedures outlined in Wang et 
al. (2004). In particular, we write the 2PCF as 



^(r, z) = fih(r, z) + £,2h{r-, z) , 



(27) 



where ^ih represents the correlation due to pairs of galax- 
ies within the same halo (the "1-halo" term), and ^2h 
describes the correlation due to galaxies that occupy dif- 
ferent halos (the "2-halo" term). 

In our model, we assume that the radial number den- 
sity distribution of the satellite galaxies, ns(r|M, z) fol- 
lows the same NFW profile as a dark matter halo of mass 
M at redshift z out to the virial radius, rvir, defined as 
the radius inside of which the average density is Avir 
times the critical density. To good approximation 



Avir — IStT 



60a; - 32a;^ 



(28) 



with X = ^^{z)-l (Bryan & Norman, 1998). Note that 
in reality, the distribution of satellite galaxies may not 
follow an exact NFW profile. However, as tested in Yang 
et al. (2004) by populating satellite galaxies according 
to the real mass distributions in dark matter halos, the 
change produced in the 2PCFs is quite small, less 20% on 
small scale in real space. For the projected 2PCFs which 
will be used in this work, the impact is even smaller, less 
tha 10%. Furthermore, our projected 2PCFs are to be 
estimated on scales larger than 0.1 /i~^Mpc, where the 
results are not sensitive to the assumption on the concen- 
tration of galaxy distribution in dark matter halos (see 



The Galaxy - Dark Matter Connection 



7 



Yang et al. 2005). Thus the foible of using NFW pro- 
file is negligible in comparison to the uncertainties in the 
observational data and in the halo model itself. For con- 
venience, in what follows we use the normalized number 
density distribution, Us{r\AI,z) = ns{r\M, z)/ {Ns\M, z), 
so that ^ 

47r / Us{r\M,z)r^ dr = 1, (29) 

where rvir = ?'vir(Af, 2) is the virial radius of the halo 
of mass M at redshift z. The 1-halo term of the galaxy 
correlation function can then be written as 

2 f°° 

6h(r,z) = =27T / nh(M,z)(A^pair|M,z) 
ng[zj Jq 

f{r\M,z)dM, (30) 

where (A'pairl , z) is the mean number of distinc- 
tive galaxy pairs in a halo of mass M at redshift z, 
f{r\M, z)47rr^Ar is the fraction of those pairs that have 
separations in the range r± Ar/2, and ng(z) is the mean 
number density of galaxies at redshift z, given by 

/•oo 

ng(z)^ ni,{M,z){N\M,z)dM . (31) 
^0 

For a given (M, z), the mean number of pairs as function 
of separation, (A^pair) /(?"), can be divided into contri- 
butions from central-satellite pairs and satellite-satellite 
pairs: 

(A^pair)/(r) = (iVcs)".(r) + {N,s)fsir) , (32) 

where fs{r) is the satellite galaxy pair distribution func- 
tion within a dark matter halo of mass M at redshift 



/s(r)=27r/ Us{s)s'^ds Us{\s + r\) sin d0 , (33) 



with |s + r| = (s^ 



2sr cos 9Y^^ . The number of 



central-satellite pairs is 

{N,,\M,z) = {N,\M,z) (iVs|M,z) 



(34) 



For the number of satellite-satellite pairs we assume that 
follows Poisson statistics so that 



{N,,\M,z) = (iVs|M,^) 



(35) 



The 2-halo term of the 2PCF for galaxies can be writ- 
ten as, 



^2h{r, z) = [/oxcl?-, z)] ^dm{r, z) 



(36) 



where ^.dmir, z) is the non-linear 2PCF of dark mat- 
ter particles at redshift z, which we obtain by Fourier 
Transforming the non-linear power spectrum of Smith ct 
al. (2003), and 



fe^c{r,z): 



1 



fH,{M,z){N\M, z) 



ng{z) JO 
bi,{M,z)U{r\M,z)dM . 



(37) 



Here we have taken into account the halo-halo exclusion 
effect, in that the 2-halo galaxy pair can not have an 
average distance smaller than rc->cc{M, z): 

[/(r|M,z) = |? L:'^'^'^^^'^''^ (38) 



As shown in Wang et al. (2004), this method of comput- 
ing the 2PCFs is accurate at the few percent level as long 
as rcxc{M^ z) = 2rvir(M, z), which is the value we adopt 
throughout. 

Observationally, the direct measurement is not the 
real-space 2PCF because of redshift space distortions. 
Instead one measures the projected 2PCF, Wp(rp), which 
is related to the real space 2PCF as 



Wp{rp) = 2 ^(^rl+rl) dr^ , 



(39) 



where the comoving distance has been decomposed into 
components parallel (rjr) and perpendicular (rp) to the 
line of sight. In practice, the integration is only per- 
formed over a finite range of Vt^. Following Wang 
et al. (2007), we use a maximum integration limit of 
'"TT.max = 40 /i^^Mpc (see Appendix A). 

3.3. Redshift Dependence 

The formalism described above allows us to model the 
stellar mass functions and correlation function of galax- 
ies once a set of parameters is adopted to specify the 
CSMF, ^ciM^\M,z), and the dynamical evolution func- 
tion, P{m^,\m^.^a, Za',ma', Ma^rj, z). At any particular red- 
shift, our model for these quantities is fully described by 
the following seven free parameters: qj -^Ij ctj Pi 
Pt and c. In order to describe the evolution of the galaxy 
distribution over cosmic time, we need to specify how 
each of these parameters changes with redshift. 

Throughout this paper we assume the following red- 
shift dependence for our model parameters: 

log[M,,o(^)] =log(A/*,o) + liz 
log[Mi(z)]=log(A/i) -f^72Z 
a{z)^a + 73Z 
log[/3(z)] = min[log(/3) + 74Z + 75^', 2] (40) 
f7c(z) = max[0.173,0.2z] 
Pt{z)=pt 
c(z) = c . 

These functional forms have been obtained using trial- 
and-error whereby we tried to minimize the number of 
parameters while achieving a good fit to the overall data 
(see Section U below) . As we have tested and as pointed 
out in Moster et al. (2010), adz) cannot be well con- 
strained by the SMFs alone. The form adopted above for 
<Jc{z) is based on the amount of true scatter at low red- 
shift obtained by YQ9 and the uncertainty in the stellar 
mass measurements at high redshifts (e.g., D05; see also 
Conroy et al. 2009). Note that the redshift dependence 
is relevant only at high redshift when cTc is dominated by 
the uncertainty in the stellar mass measurements, and so 
we use o'c(z), instead of (Tc(za), for CTc in Eq. [551 The 
model is thus specified by a total of 11 free parameters, 
four to describe the CSMF at z = 0, five to describe their 
evolution with redshift, and two (c and pt) to describe 
the evolution of satellite galaxies. Note that we assume 
that both c and pt are independent of redshift. As we 
demonstrate below, this assumption is perfectly compat- 
ible with current data. However, it is very simplified 
and needs to be re-evaluated once better data becomes 
available. 
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4. OBSERVATIONAL CONSTRAINTS 

The goal of this paper is to constrain the CSMF, 
$c(-M*|M, z), and the dynamical evolution function, 
P(m*|TO,^a, Za', rua', Ma', rj, z), and in particular their evo- 
lution with redshift, using observational data. In this 
section we summarize the observational data to be used 
to constrain our model. All our low redshift data are 
obtained from the SDSS DR7 (Abazajian et al. 2009). 

4.1. Low- Redshift Data 
4.1.1. The Total Stellar Mass Function 

For the low redshift stellar mass function (SMF) we 
use the results that are obtained using the same method 
outhned in Y09b, but here updated to the SDSS DR7 
(Abazajian et al. 2009). Note that in these measure- 
ments the galaxy stellar masses are estimated using the 
model of Bell et al. (2003) in which a Kroupa (2001) IMF 
was adopted (Borch et al. 2006). The resulting SMFs, for 
galaxies of different colors and of different group member 
types, are provided in Appendix B for reference. Here for 
our purpose, we only use the SMF of all galaxies, which 
is shown in the upper panel of Fig. [2] as open circles with 
error bars. For comparison, the solid and dashed lines 
show the best Schechter fits for the SMFs of our measure- 
ments and the one obtained by Panter et al. (2007) from 
the SDSS DR3 (Abazajian et al. 2005), respectively. The 
open squares correspond to the SMF obtained by Li & 
White (2009) from the SDSS DR7 as well but for stellar 
masses obtained by Blanton & Roweis (2007). To have 
a clearer vision of the difference between these measure- 
ments, we show in the lower panel of Fig. [2] the ratio of 
the SMFs which are normalized by our best fit Schechter 
form SMF. There are - 50% and - 20% difference at 
the high and low mass ends between these three SMFs, 
which most likely arise from differences in (i) the meth- 
ods used to derive the stellar masses and (ii) the differ- 
ent data releases used. As pointed out in Li & White 
(2009), srarting from similar IMF, the systematic differ- 
ence of stellar masses obtained by different approaches 
are at O.ldex to O.Sdex at high to low mass ends. For 
the systematic errors, it is very difficult to come up with 
a realistic covariance model for the stellar mass function, 
as it is not clear the uncertainties in the assumed IMF 
and spectral synthesis model. In the simplest assump- 
tion, the covariance can be a constant shift in the stellar 
mass, where our results can be scaled to other IMFs and 
spectral synthesis models. Here, for simplicity, to mimic 
such kind of systematics, we assume there is at least 20% 
error on each data point of our SMF. 

4.1.2. The conditional stellar mass function 

In addition to the (global) SMF, we also use the con- 
ditional stellar mass function (CSMF) to constrain our 
model. The CSMFs have been obtained by Y09b from 
the SDSS DR4 group catalogues constructed by Yang et 
al. (2007) using the halo-based group finder of Yang et 
al. (2005a). Here we update the CSMFs using the latest 
SDSS DR7 group catalogued The halo masses provided 
in the group catalogues are obtained using the halo- 
abundance matching method (see Section 4.1 in Yang et 

* see \http://gax.shao.ac.cn/data/Group.htrnI\ 




7 8 9 10 11 12 
log [M^ / h-2Mj 

Fig. 2. — Upper panel: the stellar mass function (SMF) of galax- 
ies at low redshifts {z ~ 0.1). Open circles with errorbars are 
the results we obtained from the SDSS DR7, and the solid line 
is the corresponding best-fit Schechter form. For comparison, the 
dashed line shows the best-fit of the SMF obtained by Panter et 
al. (2007), and the squares with errorbars are the results obtained 
by Li & White (2009). Lower panel: similar to the upper panel, 
but here for the ratios between the SMFs and our best-fit Schechter 
form. 

al. 2007) where halo mass function for the WMAP7 cos- 
mology is used. The results of the CSMFs are provided 
in Appendix B for reference. 

4.1.3. The two-point correlation function 

We have measured the projected 2PCFs for SDSS DR7 
galaxies of different stellar masses, estimated using the 
model of Bell et al. (2003), i.e., these are the same stellar 
mass estimates as used for our SMF and CSMF. A de- 
tailed description of how we obtained these 2PCFs can 
be found in Appendix A. We note that similar measure- 
ments have already been carried out by Li et al. (2006) 
using data from the SDSS DR2 with stellar masses as 
given by Kauffmann et al. (2003). As we have tested, 
our results agree with those obtained by Li et al. (2006) 
at 1-CT level. Here we make use of our own measurements 
with full covariance matrixes to constrain the CSMFs. 

4.2. High-Redshift Data 

In addition to the low-redshift data discussed above, 
we also use data at high redshifts. Numerous studies 
have measured SMFs at high redshifts using different 
data sets. Unfortunately, differences among the SMFs 
obtained by different authors are still quite large. In this 
paper we therefore use two different sets of data, and 
use the differences in the resulting models as a method 
to gauge the model uncertainties arising from the uncer- 
tainties in the data. 

The first data set that we use are the SMFs of Perez- 
Gonzalez et al. (2008, hereafter PG08) obtained from 
Spitzer data for galaxies in the redshift range from 
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Fig. 3. — The stellar mass functions of galaxies in different redsfiift bins as indicated in eacfi panel. The filled circles with error-bars are 
results obtained from the Spitzer observation by Perez-Gonzalez et al. (2008, PG08). The open squares with error-bars are results obtained 
in this work from the SDSS DR7. The filled triangles with error-bars are results obtained by Drory et al. (2005, DOS), for which we have 
combined their data for FORS Deep and the GOODS/CDF-S Fields. The redshift range labeled in the upper right-hand corner of each 
panel is for the results of PG08, while the range labeled in the lower left-hand corner is for the results obtained from SDSS DR7 (first 
panel) and DOS (other panels). Here the stellar masses of PG08 and DOS are divided by a factor of 1.7 to take account the systematic 
differences of using different IMFs. 



^ z < A. Note that in their study, PG08 assumed a 
Salpeter (1955) IMF which is different from the Kroupa 
IMF used in the local universe for the SDSS stellar 
masses. As suggested in PG08, the stellar masses based 
on the Salpeter (1955) IMF are systematically larger by 
a factor of ^ 1.7 than that based on the Kroupa IMF. 
For consistency, we therefore devide all the stellar masses 
provided in PG08 by this factor. The resulting data are 
shown as solid dots with errorbars in Fig. [3l For com- 
parison, the open squares in the upper left-hand panel 
correspond to the local SMF obtained from the SDSS 
DR7 in this study. Note that the latter is in general con- 
sistent with that of PG08 except at low mass end, which 
we believe to be a result of cosmic variance and/or the 
stellar mass estimator. In what follows we therefore ig- 
nore the SMF of PG08 for the < z < 0.2 redshift bin. 

The second data set for high redshift SMFs is that 
of Drory et al. (2005; hereafter D05), which is based 
on deep multicolor data in the FORS Deep Field 
and the Great Observatories Origins Deep Survey- 
South/Chandra Deep Field-South (GOODS/CDF-S) re- 
gion. For our investigation, we combine the two separate 
measurements of these two regions, taking their average 
and using the maximum of their original error and the 
difference between the two as the new error. For the 



few data points for which only FORS Deep Field data 
is available, we update their errors using the percentage 
errors of the nearest data points to partly take account 
of systematic errors. Again, the stellar masses obtained 
in D05 assumed a Salpeter (1955) IMF, and we divide 
the stellar masses by a factor of ~ 1.7 to correct them to 
the Kroupa IMF. The resulting data are shown in Fig. [3] 
as filled triangles with errorbars. Note that there is large 
discrepancy between the PG08 and the D05 results, es- 
pecially at high redshift. Since it is difficult to decide 
which result is more reliable, we will perform two sepa- 
rate analyses using either of these data sets. Finally, as 
we did for the low-redshift data, we also assume a min- 
imum of 20% error for each data point of the SMFs to 
account for other possible uncertainties in the data. 

Two-point correlation functions for galaxies of different 
stellar mass bins are available at intermediate redshifts 
(e.g., Meneux et al. 2008 using VVDS; Meneux et al. 2009 
using zCOSMOS). However, as pointed out in Meneux et 
al. (2009), because of the small sky coverage of the sur- 
veys, these measurements suffer significantly from cosmic 
variance. We therefore decided not to use these measure- 
ments to constrain our model. Instead we shall only use 
them for comparison with our model predictions. 
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TABLE 2 

Model Abbreviations 



Model 


Data constraint 


c 


Pt 


FITO 


SMF 







FITl 


SMF 


1 


oo 


FIT-2PCF 


SMF+2PCF 


free 


free 


FIT-CSMF 


SMF+CSMF 


free 


free 


SMFl 
SMF2 


SMF=SDSS DR7+PG08 
SMF=SDSS DR7+D05 







Note. — SDSS DR7: the stellar mass function at z = 0.1 wo 
obtained from SDSS DR7 using the method of Yang et al. (2009b); 
PG08: the stellar mass funetions obtained by Perez-Gonzalez et 
al. (2008) from Spitzer, excluding the data at z — 0.1; DOS: the 
combined stellar mass functions obtained by Drory et al. (2005) 
using deep multicolor data in the FORS Deep Field and the 
GOODS/CDF-S region. 



4.3. The Fitting Method 

To constrain model parameters using observational 
data, we use a Monte-Carlo Markov Chain (hereafter 
MCMC) to explore the likelihood function in the multi- 
dimensional parameter space (see Yan, Madgwick & 
White 2003 and van den Bosch et al. 2005b for simi- 
lar analyses). In order to check whether different data 
sets are mutually compatible, we carry out two different 
analyses. The first one, referred to as FIT-2PCF, uses 
the SMFs together with the 2PCFs [wp{rp)] of galaxies 
at low-z. The second, referred to as FIT-CSMF, uses the 
SMFs at different redshifts together with the CSMFs at 
low z as constraints. The corresponding are defined 
as 

x1pcf=x'(<J^)+x'K), (41) 



and 



XCSMF = + X^(*CSMf) , 



(42) 



respectively. Here x^('I'), X^i^p), x^('J'csmf) are given 

by 



i=l j 



$mod(M*j, ^i) - ^ohsi^Uj, Zi) 



1 2 



a.s,{M^^j,Zi) 



and 



with 



(43) 



X\wp) A«;pj(C-i),- fcA«;p,k , (44) 

i=l j,k 



AWpj = Wp,mod(?'pj|Af*,i) - Wp,obs(?'pj|Af»,i) , (45) 



and 



CSMFj 

Nc 



1=1 J 



^mod(M,,j-|M,) - ^oUM*,j\Mi 

crcSMF{M,j\Mi) 



(46) 



and 



In the above expressions, (T$(M, j, z^), 
acSMF{M^,j\Mi) are, respectively, the observational 
errors for the SMF and CSMF, while (C"^)j,fc is the in- 
verse of the covariance matrix of the observed projected 
2PCF. The integers N^, and Nc are the number of 
SMFs (for different redshifts), the number of projected 



2PCFs (for different stellar mass bins), and the number 
of CSMFs (for different halo mass bins), respectively. 
The factor 1/Nc in Eq. (|46)) is somewhat arbitrary, 
which is included to reduce the weight of the CSMFs, as 
both the CSMFs and the low-z SMF are based on the 
same data and thus not completely independent. 

In addition to FIT-2PCF and FIT-CSMF, in which 
all eleven parameters are included in the fitting, and 
both the SMFs and 2PCFs (or CSMFs) are used as con- 
straints, we also consider two extreme cases where pt is 
assumed to be either or oo. The first case (referred to 
as FITO in the following) with pt — effectively assumes 
that all satellite galaxies are stripped as soon as they are 
accreted (i.e., there are no satellite galaxies). The second 
case (referred to as FITl) assumes pt = oo and c = 1, so 
that a satellite is never disrupted and evolves in such a 
way that, at the redshift z in question, it has the same 
stellar mass as a central galaxy in a host halo with mass 
equal to the halo mass of the satellite at its accretion 
time. In these two extreme cases we only use the SMFs 
at different redshifts as constraints. These extreme cases 
are used to gauge the importance of the satellite popu- 
lation in our modeling (see also Yang et al. 2009a; Wang 
& Jing 2010). 

As mentioned above, there are large differences be- 
tween the SMFs obtained by PG08 and D05 at high 
redshifts, and we therefore will use them separately to 
constrain our model. More specifically, we consider two 
sets of models. The first uses the high-z SMFs of PG08 
and is refered to as SMFl; the second uses the high-z 
SMFs of D05 and referred to as SMF2. For reference. 
Table [5] lists all our different models, including the ab- 
breviations that we use to refer to them. 

We start our MCMC from an initial guess and allow a 
'burn-in' of 1,000 steps for the chain to equilibrate in the 
likelihood space. At any point in the chain we generate a 
new set of model parameters by drawing the shifts in the 
eleven free parameters from eleven independent Gaussian 
distributions. The Gaussian widths are tuned so that the 
average accepting rate for the new trial model is about 
0.25. The probability of accepting the trial model is 



accept 



1.0 

exp[-(x,i 



if X 

xL)/2] ifx'" 



2 < 

new Aold 
new — Aold 



(47) 

with the x^ defined in Eq. (|4T|) or in Eq. (|42|) . 

We construct a MCMC of 500,000 steps after 'burn- in', 
with an average acceptance rate of ~ 25 percent. Before 
the further treatment of the MCMC, we have tested its 
convergence using the 'convergence ratio' r as defined in 
Dunkley et al. (2005). In all cases in consideration, we 
observe r < 0.01 is achieved for each parameter. In order 
to suppress the correlation between neighboring models 
in the chain, we thin the chain by a factor 50. This re- 
sults in a final MCMC consisting of 10, 000 independent 
models that sample the posterior distribution. As an il- 
lustration, we show in Fig. |4]the fiowchart of our filing 
process. Based on these 10, 000 independent models, we 
obtain the best parameter set that has the smallest x^- 
The 68% confidence level (CL) around this set of model is 
obtained by simply sorting the models according to their 
X^, and finding models that have values smaller than 
that corresponding to the 68% CL. As an illustration, 
the outer to inner contours in Fig. [5] plot the projected 
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The Conditional Stellar Mass Function: <1>(/W* \M, z) = ^c^W* \M , z) + <t>s(m* \M , z) 



c!>s(m, |/W, z) = f{<^c(m* | , Zj), n^„i.,(m3, z^l/W, z), pt, c) 



Stellar Mass Functions 
Projected 2PCFs 

M/p{''p , Z = 0.1) 

CSMFs of centrals 
(t>c(M* |/W, z = 0.1) 

CSMFs of satellites 
*s(m*|/W,z = 0.1) 




Halo mass function 
nh(/W,z) 



Halo density profile 



Halo bias model 
bi,(/W,z) 



Subhalo accretion model 
"subtnis, zaiM, z) 



Run MCMC to constrain the parameters: {M^,^0; 0^ 7i 1 1'2: 73 > 74? 75 i Pti '^l 



Fig. 4. — A flowchart showing how we constrain the paramenters of CSMFs. The fllled bars on the left-hand side hst the observational 
quantities used as constraints, while those on the right-hand side list the ingredients to model these quantities. We use MCMC to obtain 
both the best fit parameters and their confidence ranges. The MCMC analysis is performed separately for both FIT-2PCF and FIT-CSMF. 



2-D boundaries in the parameter space that enclose 95%, 
68%, 30%, 10% and 1% of the models with the smallest 
values, with the best-fit values indicated by a cross. 
The marginalized 1-D distribution of each parameter ob- 
tained from all the 10, 000 models is also shown in the 
plot. As one can see, strong covariance exists in some 
parameter pairs, such as /3 and Mi, 71 and 73, and 74 
and 75. Overall, most parameters are well constrained, 
with their mean values much larger than the correspond- 
ing dispersions. For reference we list in Tables |3] and |4] 
the best-fit parameters together with their 68% CLs. IDs 
1 to 4 correspond to FITO, FITl, FIT-2PCF and FIT- 
CSMF, respectively, all assuming WMAP7 cosmology. 
Results are listed separately for cases using SMFl and 
SMF2. For the other four cosmological models listed in 
Table [T] we only give the best-fit parameters and their 
68% CLs for FIT-CSMF and FIT-2PCF models (IDs 5- 
12). The details of all these models are discussed in the 
following section. 

5. RESULTS 

The upper left-hand panels of Figs. [6] and [7] show the 
fits to the observed SMFs at low-z using SMFl and 
SMF2, respectively. In both cases we adopted the fiducial 
WMAP7 cosmology. The long-dashed and dotted lines 
show the results for the two extreme models, FITO and 
FITl, respectively. Clearly, even these two extreme mod- 
els can accurately match the observed SMFs (as judged 
by their reduced values), indicating that the SMFs 
alone cannot constrain the values of pt and c. This is ex- 



pected from the fact that, as mentioned above, the SMFs 
are dominated by central galaxies. Note that for SMF2, 
the x^{^) values are significantly larger than those for 
SMFl. The main reason for this is the much smaller er- 
ror bars in SMF2 and the large difference between model 
and data at the low mass end for the 0.25 < z < 0.75 
redshift bin (see the second panel of Fig.[71). 

The best fit parameters of the FITO and FITl models 
are used to make predictions for the projected two-point 
correlation functions, Wp(rp), in different stellar mass 
bins, and the results are shown in Fig. |8] (for SMFl) 
and Fig. |9] (for SMF2), again as long-dashed and dot- 
ted lines, respectively. Clearly, the model predictions 
for FITO and FITl are very different, especially on small 
scales where the correlation function is dominated by the 
1-halo term. Since FITO and FITl represent the mini- 
mum and maximum contributions of satellite galaxies, 
respectively, the correlation functions predicted by these 
two models should bracket the predictions of more realis- 
tic models. The large difference between FITO and FITl 
clearly demonstrates that the observed Wp{r-p)^ especially 
for < 1 /i^'^Mpc, can provide useful constraints on 
the satellite population, which in our model is character- 
ized by the parameters pt and c. Note that the observed 
u;p(rp) falls nicely between the predictions for our two 
extreme models, suggesting that a more realistic treat- 
ment of satellite galaxies should be able to fit the data. 

When separated into central and satellite components, 
the CSMF, $(M*|M), provides direct information re- 



TABLE 3 

The best fit parameters of the CSMF based on SMFl. 



ID 


cosmology 


log A/,,0 


log A/i 


a 


/3 


71 


72 


73 


74 


75 


Pf 


c 




FIT-type 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 


(14) 


(15) 


1 
2 

3 

4 


WMAP7 


io.2i«;?;; 

10.12«;i; 
in 2's+° °^ 

iU.^a_Q Qg 

in iq+0.08 


10.73t«:i^ 

io.821e;:16 

70+0.06 

iU./ci_(j Q5 


n qq + 0.02 

"^■■-'■^-0.02 
n 00+0.02 

'J-^y-0.02 


5 K4+1.95 
^■'^^-i.io 
9.401^-66 

IQ+0.49 

^■J^y-0.59 
^■15-2.33 


-0.67l^J2 
-n 'ii+o is 

'^■"^-0.19 


-n 27+°-''^ 

'-'■^'-0.06 

r, 97 + 0.06 

'-'■^'-0.06 

-n 2fi+° °3 

'-'■^"-0.05 

-n 24+°-''5 


'J-J^^-0.08 
^•J^°-0.07 

n 10+0.07 


0.14lSJ-^ 

-o.o3±s;:i^ 
-o.o7is;:i^ 


'J-^3-0.08 

n 90+0. 08 

n qO+0.09 

'^••^■J-0.06 

n .„+o.io 

'-'••^■'-0.09 




00 

, 01 +0.23 
^^■■'^^-0.15 

n 09+0.10 

'J-°^-0.08 


1.00 

'^•"^-0.58 
r\ '7Q + 0.23 
^- '■J-0.24 


35.10 
38.91 
113.33 
102.09 


FITO 

FITl 
FIT-2PCF 
FIT-CSMF 



(SO. 



5 


WMAPl 


10, 


9,,+0.04 
■^•J-0.07 


10. 


0.+0.02 

■84_o.i4 


0, 


90+0.01 
■^^-0.01 


7. 


1 1 +1.64 

-'--'-0.25 


-0, 


,62«: 


16 

,16 


-0, 


90+0.05 

■^•'-0.04 


0. 


10+0. 06 
-'"-0.06 


_o 19+0.10 
'^•-'^-0.12 


0. 


qt; + 0.07 
■''•-'-0.07 


1, 


-19+0.16 
■-'^-0.29 


0, 


47+0.53 

■^' -0.23 


107.07 


FIT-2PCF 


(D 


6 




10, 


99+0.03 
■^^-0.07 


10. 


70+0.03 

■'^'-o.io 


0, 


90+0.01 
■^"-0.01 


7. 


gg+1.51 
'^°-0.11 


-0, 


■59t^: 


14 
,16 


-0, 


99+0.06 
■^^-0.05 


0. 


-,7+0.07 
^'-0.05 


-o.i4t«:J^ 


0. 


.,0+0.06 

■•'"-o.io 


0, 


05+0.10 
■''''-0.10 


0. 


9O+0.19 
■^^-0.15 


103.77 


FIT-CSMF 




7 


WMAP3 


10, 


97+0.06 
■^'-0.03 


10. 


79+0.09 
■'^-0.09 


0, 


90+0.01 
■^"-0.02 


7. 


qQ + 1.20 
'-'^-0.81 


-0, 


.63«: 


18 
,13 


-0, 


9Q+0.05 
■^^-0.04 


0. 


-10+0.05 
^"-0.07 




0. 


97+0.07 
■^'-0.06 


1, 


97+0.44 
■^' -0.20 


0, 


04+0.06 
■^^-0.38 


151.37 


FIT-2PCF 




8 




10, 


91+0.04 
■^-^-0.04 


10. 


■"^-0.05 


0, 


9Q+0.01 
■^"-0.01 


7. 


gg + 0.59 

'^'J-0.88 


-0, 


,62«: 


,20 
,15 


-0, 


97+0.04 
■^'-0.04 


0. 


1Q+0.04 
^^^-0.08 


10+"-" 

'J-^'J-0.08 


0. 


9c;+0.06 
■^■'-0.07 


0, 


77+0.08 
■"-0.06 


0. 


Q7+0.03 
■^'-0.13 


100.64 


FIT-CSMF 




9 


WMAP5 


10, 


99+0.05 
■^^-0.06 


10. 


7n+0.06 
■"^-0.09 


0, 


9Q+0.02 

■^'^-o.oi 


8. 


10 + 1. 04 
-"'-0.68 


-0, 


.59^0; 


14 

,18 


-0, 


9tr+0.04 
■^■^-0.05 


0. 


-,0+0.06 
^"-0.06 


-0 01+''-" 


0. 


9Q+0.07 
■^^-0.07 


1, 


10+0.20 
■^"-0.12 


0, 


QK + 0.05 

■y^-0.51 


120.92 


FIT-2PCF 




10 




10, 


99+0.03 
■^^-0.06 


10. 


70+0.05 
■'■'-0.07 


0, 


9Q+0.02 
■^^-0.01 


7. 


nK+0.89 
■'^^-0.57 


-0, 


.eot": 


,20 
,16 


-0, 


9C+0.05 
■^"-0.04 


0. 


-,0+0.06 
^"-0.08 


0.03+n? 


0. 


90+0.07 
■^"-0.06 


0, 


00+0. 08 
■"'-'-0.08 


0. 


go+0.18 
,DU_Q 21 


106.20 


FIT-CSMF 




11 


Millennium 


10, 


97+0.04 
■^'-0.05 


10. 


74+0.08 
■'^-0.10 


0, 


97+0.02 
■^'-0.01 


7. 


ig+1.31 
-'"-0.80 


-0, 


.67«: 


22 
,11 


-0, 


90+0.06 
■^■'-0.05 


0. 


9n+0.04 
^^^-0.08 


_OQ9+0.13 

'J-'-'^-o.io 


0. 


0-, +0.06 
■''-'-0.08 


0, 


04+0.43 
■^^-0.07 


0, 


Q7+0.03 
■^'-0.74 


98.43 


FIT-2PCF 




12 




10, 


99+0.05 
■^^-0.06 


10. 


7n+o.io 

■'"-0.08 


0, 


90+0.02 
■^"-0.01 


7. 


9C+O.94 
■^^-1.14 


-0, 


■57t°: 


16 

,18 


-0, 


9-1+0.04 

■^-'-o.oe 


0. 


-lfi+0.06 
^"-0.06 


QQ,+0.13 

"•'-'^-0.08 


0. 


.10+0.05 
■''"-0.07 


1, 


QQ+0.10 

■"'J-0.08 


0, 


qtr + 0.16 
■'-'''-0.13 


98.28 


FIT-CSMF 





Note. — Column (1): model ID. Column (2): cosmology used for fitting. Column (3-13): the best fit parameters for the CSMF. Note that M*^o is in units of h "^Mq and M\ in units of 
h~^M(7), and log is the 10 based logarithm. Column (14): the defined in Eq.l l41| i or ( 142D . Column (15): the type of FIT used to constrain the models. Here results are obtained by using 
SMFl . 
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Fig. 5. — The best fit (cross or vertical fine) and the projected distribution of the parameters in 2-D or 1-D space. The outer to inner 
contours in the 2-D plane correspond to the boundaries that enclose 95%, 68%, 30%, 10% and 1% models with the smallest values, 
respectively. The 1-D distributions are the marginalized distributions of individual parameters obtained from all 10, 000 models. 



TABLE 4 

The best fit parameters for the CSMFs based on SMF2. 



ID 


cosmology 


log 


logA/i 


a 


/3 


71 


72 


73 


74 


75 


Pt 


c 




FIT-type 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 


(14) 


(15) 


1 
2 
3 
4 


WMAP7 


10.45«;?I 
10.45+0-;? 

in sfi+0 05 


11 17+0-11 
11-1'-0.19 
11 43+0.15 

11 21+0-" 
11-^1-0.09 

ii.neini 


n 25+0 03 
n 97+0.03 

^-^'-0.02 

n 25+0 01 

"■^■^-0.03 

n 27+0 01 


q r,q+0.78 

■^-^•J-0.37 
q oc:+0.92 

4 nn+0-54 

4 34+0-96 
*-'^*-0.52 


-n.98tO-22 
-n.79t0-23 

-n Q3+0-14 

-n qfi+O-13 


-n 25+0-07 

'-'-^''-0.05 

-n 24+0-10 
-n 26+0-07 

'-1-^0-0.03 

r, r,q + 0.05 

'-'-^'-'-0.06 


0.42+0-0? 
0.3810-0 

SQ+O-ll 
^•''^-0.05 
41+0-07 
^•^1-0.08 


-0 n7+o-i3 

^'■'J' -0.09 

-0 16+0-09 
^•1°-0.09 
-0 15+0-07 

-0 11+0-11 
'J-11-0.08 


ni+0-05 

'^•'-'1-0.01 

02+0-04 
'^•*J^-0.02 

04+0-02 

ni+0-05 

^•'Jl-0.07 




oo 

, 1O+0.21 
l-lS-0.23 

0.88^0-07 


1.00 

0.911°:°^ 

0.9810:02 


56.03 
52.91 
128.54 
131.21 


FITO 

FITl 
FIT-2PCF 
FIT-CSMF 



CD 



5 


WMAPl 


10, 


46+0-06 


o(,+0.09 


0, 


nc+0.02 
■^3-0.02 


3 6Q+0-85 


-0, 


07+0.10 

■°'-0.18 


-0, 


9n+0.06 

■^^-0.06 


,-, q7+0.09 
'J-'^' -0.05 


n -,0+0.08 


04+0-03 


0, 


Q7+0.23 
■^'-0.14 


0.94l°:°« 


123.92 


FIT-2PCF 


6 




10, 


qQ+0.04 
■•^^-0.06 


-|i -17+0.09 

ii-i'-o.io 


0, 


nf;+0.02 
■^"-0.01 


4 34+0.62 
*-'^*-0.54 


-0, 


Qq+0.16 
■^■^-0.12 


-0, 


9-, +0.05 
■^1-0.09 


n qo+0.06 


-0 17+0-10 
'J-l'-0.08 


ns+o-os 

l'-'J''-0.04 


1, 


nq+0.13 
■'-'■^-0.09 




131.82 


FIT-CSMF 


7 


WMAP3 


10, 


49+0.03 
■^^-0.05 


11 16+0-05 

ii.iD_o gg 


0, 


94+0.02 
■^*-0.02 


3 0-1+0.45 
■^■^1-0. 37 


-0, 


Qq+0.13 
■^■^-0.14 


-0, 


90+0. 07 
■^"-0.04 


41+0-07 
'^•^1-0.06 


-0 11+0-10 

'J-11-0.08 


n3+o-02 

'J-'-' -^-0.08 


1, 


qn+0.25 
■■^'-'-0.17 


n qq+0-01 


167.19 


FIT-2PCF 


8 




10, 


,7+0.04 
■■^'-0.06 


10.94lO:J| 


0, 


9f.+0.02 
■^"-0.01 


4 09+0. 61 


-1, 


nq+0.19 
■'-'■^-0.10 


-0, 


9C+0.06 

■^"-0.04 


n 43+0-06 


-0 04+0-11 


00+0-03 

U.UU_o 07 


0, 


otr+O.lO 

■"■^-o.oe 


i-ool«:~ 


126.71 


FIT-CSMF 


9 


WMAP5 


10, 


cfi+0.03 
■"^"-0.14 


11.3110-0^ 


0, 


99+0.05 
■^^-0.01 


3 54+0-75 
"^■"^^-0.33 


-0, 


qn+O-iO 


-0, 


9tr+0.06 
■^°-0.05 


42+0-07 
"•^^-0.07 


-0 16+0-10 
'J-lO-0.08 


n3+o-03 

U.UCi_o 06 


1, 


9C+0.21 

■^3-0.17 


0.98i0:07^ 


135.77 


FIT-2PCF 


10 




10, 


qc+0.05 
■■^"-0.06 


ii.oetni 


0, 


97+0.01 
■^'-0.01 


4 34+0.96 
*-'^*-0.52 


-0, 


qfl+0.13 
■^"-0.19 


-0, 


9q+0.05 
■^■^-0.06 


41+0-07 
'J-^1-0.08 


-0 11+0-11 

'J-11-0.08 

-0 in+0-05 
'J-i'J-o.io 


ni+0-05 

'^•*Jl-0.07 


0, 


00+0. 07 
■O'^-O.IO 
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The same as Table[3l but here SMF2 is used instead of SMFl. 
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Fig. 6. — The predicted stellar mass functions (SMFs) of galaxies in different redshift bins compared with the observed SMF of SDSS 
DR7 and those of PG08 (other panels). In each panel, the long-dashed, dotted, dashed and solid lines, which are quite similar, show the 
predictions of models FITO, FITl, FIT-2PCF and FIT-CSMF, respectively, for the case of SMFl and assuming the WMAP7 cosmology. 
In each panel, the dot-long-dashed lines are the predictions for the "no-evolution model" in which we have assumed that the central galaxy 
- halo relation does not evolve and is the same as that at 2 = 0.1. Here results for (7c{z) = 0.173 (smaller at massive end and at high 
redshift) and crc(z) = max[0.173, 0.2z] are plotted together for comparison. 




log [M^ / h-2Mg] log [M, / h-2Mg] log [M^ / h-2Mg] log [M^ / h-^M^] 

Fig. 7. — Similar to Fig. |6] but here SMF2 has been used as constraints. The data shown (symbols) are the stellar mass functions of 
SDSS DR7 (first panel) and DOS (other panels). 
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Fig. 8. — The projected 2PCFs of galaxies in different stellar mass bins as indicated by the log[A/*/(/i ^ M0)]-vaIues in brackets. Open 
circles with error bars are the observational measurements from the SDSS DR7 described in Appendix A. The solid, dashed, dotted and 
long dashed lines show the predictions for models FIT-CSMF, FIT-2PCF, FITl and FITO, respectively, where SMFl has been used as 
constraints for the stellar mass functions. 
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Fig. 9. — Same as Fig.[8l except that here SMF2 has been used as constraints for the stellar mass functions. 
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Fig. 10. — The conditional stellar mass functions for central and satellite galaxies in halos of different masses as indicated in each panel 
(the unit of Mjj used is /i^^M©). The filled dots and open circles are the measurements from the SDSS DR7 group catalogue for central 
and satellite galaxies (see Appendix B), respectively. The error bars represent scatter among different bootstrap samples. The long-dashed, 
dotted, dashed and solid lines in each panel show the predictions for models FITO, FlTl, F1T-2PCF and FIT-CSMF, respectively. Note 
that the satellite component is absent in FITO. Here we have adopted the WMAP7 cosmology and used SMFl as constraints for the stellar 
mass functions. 



18 



Yang et al. 




8 9 10 11 12 8 9 10 11 12 



log [M^ / h-'^U^] log [M^ / h-^Mg,] 



a 

O 

i-i 
60 



00 
O 



o 



- 



-1 



T — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — p 

13.0< log Mj, ^13.3 J 




J I I I I I I I L 



9 10 11 12 
log [M^ / h-^M^] 



T — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — r 



12.K log ^12.4 




_I I I I L 



9 10 11 12 
log [M^ / h-^M^] 



Fig. 11. — Same as Fig. 1101 except that here SMF2 has been used as constraints for the stellar mass functions. 



garding the satellite fractions, and so it is also interest- 
ing to look at the model predictions for $(M, |M). The 
long-dashed and dotted lines in Fig. [10] (for SMFl) and 
Fig. [n] (for SMF2) show the z = 0.1 predictions of mod- 
els FITO and FITl, respectively. Note that by definition 
there are no satellite galaxies in FITO, so that there are 
no long-dashed lines present for the satellite components. 
For comparison, the solid dots and open circles are ob- 
servational measurements obtained from SDSS DR7 (see 
Section HTg). Note how FITO and FITl predict CSMFs 
for centrals that are very similar, and in good agreement 
with the data. However, model FITl clearly, and sig- 
nificantly, over-predicts the CSMFs of satellites, while 
FITO also clearly fails to match the data since it pre- 
dicts zero satellites. Once again, the observed CSMFs 
of satellites falls in between the predictions of these two 
extreme models, suggesting that the actual CSMF data 
can be used to constrain model parameters, especially pt 
and c. 

Finally, Fig. [1^ shows the z = 0.1 relation between 
halo mass, M, and the stellar mass of the central galaxy. 



M*_c (upper panels), and the ratio M^^c/M (lower pan- 
els). Here results for SMFl and SMF2 are shown in 
the left- and right-hand panels, respectively, while long- 
dashed and dotted lines once again correspond to the 
extreme models FITO and FITl. Note how both models 
are virtually indistinguishable at the high mass end, and 
in good agreement with the results of Yang et al. (2008) 
obtained directly from the SDSS group catalog (open cir- 
cles with errorbars), here updated according to the latest 
SDSS DR7. At the low mass end, however, FITO predicts 
M^^c that are significantly higher than FITl. The rea- 
son is simple: having no satellites forces model FITO to 
assign a larger number of massive galaxies, which other- 
wise would have been assigned to massive halos as satel- 
lites, to relatively low-mass halos as centrals. Thus a self- 
consistent modelling that takes into account the subhalo 
accretion times and the evolution of satellite galaxies af- 
ter accretion mainly impacts the Mh<_c - M relation at the 
low- mass end. Such behavior has also been noticed, e.g., 
in Moster et al. (2010) and Behroozi et al. (2010). Note 
also that the Af* c - M relations obtained from SMFl and 
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Fig. 12. — The M,_c - M relations (upper panels) and the M-t^c/M ratios (lower panels) as function of halo mass. Results are shown 
separately for SMFl (left panels) and SMF2 (right panels), assuming WMAP7 cosmology. Here the results are for galaxies at redshift 
^ = 0.1. The long-dashed, dotted, dashed and solid lines are the predictions for models FITO, FlTl, F1T-2PCF and FIT-CSMF, respectively. 
The open circles with error-bars are the direct measurements obtained from the SDSS group catalogue (e.g., Yang et al. 2008), here updated 
to the latest DR7. As an illustration, the horizontal dotted lines in the lower panels indicate the mass ratios corresponding to 100%, 10% 
and 1% of the universal mass ratio between baryonic and total matter; fib/^^m = 0.167. 



SMF2 are quite similar over the entire halo mass range. 

Next we focus on the results obtained from FIT-2PCF 
and FIT-CSMF, in which all 11 parameters are kept free, 
and the observed SMFs are combined either with the 
2PCFs (in FIT-2PCF) or the CSMFs (in FIT-CSMF) to 
constrain the model parameters. The fitting results are 
shown in Fig. [6l through [T2l as the dashed and solid lines, 
respectively. The best fit values for each of the 11 model 
parameters are listed under IDs 3 and 4 in Tables [3] (for 
SMFl) and|l(for SMF2). There are a number of points 
worth noting: 

• For a given set of SMFs (SMFl or SMF2), the re- 
sults obtained from FIT-2PCF and FIT-CSMF are 
quite similar, showing that the 2PCFs and CSMFs 
constrain the model in a similar way. 

• The best-fit values forpt fall between about 0.8 and 



1.3, which is roughly consistent with the canon- 
ical value, pt = 1.0, expected from the dynami- 
cal friction time of subhalos in numerical simula- 
tions. The value of parameter c, which is much 
larger than in most cases, indicates that there 
is evolution in the mass of satellite galaxies. The 
mass growth for low-mass satellites after their ac- 
cretion into the host halo is particularly significant, 
as indicated by the redshift dependent Ah^c ~ 
relations shown in Fig. [131 For a given set of 
SMFs (SMFl or SMF2), FIT-2PCF yields larger 
(or c) (by - 30%) than FIT-CSMF, because FIT- 
2PCF slightly over-predicts the CSMFs, as shown 
in Figs. [TO] and [H 

• Using SMFl or SMF2 leads to very similar predic- 
tions for the CSMF of central galaxies, but SMF2 
yields larger satellite fractions at the low mass end 
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(Figs, [inland [ni)- 

• Models FIT-2PCF and FIT-CSMF predict M*,^ - 
M and M^,_c/M - M relations at low-z that are 
very similar, and that lie between the predictions of 
FITO and FITl (see FigHS]). They are also in good 
agreement with the results obtained directly from 
the SDSS DR7 galaxy group catalogs. These two 
models predict a maximum for the M*^c/Af ratio 

of - 0.02/1-1 ^ j^g^^Q jj^j^gg ^ ;^Qll.8 /i-1Mq, 

indicating that galaxy formation is most efhcient 
on this mass scale (at least at low-z). Note that 
a value of M^^^/M = Qmh'^ ~ 0.03 is about a 
factor of six lower than the average baryonic to 
total mass ratio ilb/^^m — 0.167 (e.g., Komatsu 
et al. 2011) of the Universe (indicated by the up- 
per horizontal dotted lines in the lower panels of 
Fig. [T2|) . This indicates that the overall star for- 
mation efficiency is low. Furthermore, since atomic 
and/or molecular gas typically contribute only a 
small fraction of the total baryonic mass of (mas- 
sive) galaxies, this furthermore underscores that 
the process of galaxy formation has either pre- 
vented most of the baryons from cooling, or has 
somehow managed to eject them from the galaxy. 

• The central stellar mass increases with halo mass 
roughly as M" '^ at the massive end and as 

for SMFl (M* for SMF2) at the low-mass end. In 
low-mass halos, the use of SMF2 leads to more mas- 
sive central galaxies than the use of SMFl (Fig fT^ . 
This arises because the stellar mass functions in 
SMF2 are steeper at the low-mass end. 

• Although the parameters in adz) are not treated as 
free parameter in our model constraints, its impact 
is quite significant, especially at high redshift where 
the scatter is large. This is demonstrated by the 
predictions of the 'no-evolution' models shown in 
Figs. [6] and [71 Here the M*_c - M relations at high- 2 
are assumed to be the same as that at z = 0.1. Two 
cases are considered, one assuming <7c{z) = 0.173 
and the other adz) = max[0.173, 0.2z]. Note that 
these two assumptions give very different predic- 
tions for the high-z SMFs at the massive end, 
demonstrating that adz) needs to be treated prop- 
erly in order to obtain reliable predictions for the 
high-z SMFs at the massive end. 

• Finally, as demonstrated by the similar values 
for different ACDM models, we note that the clus- 
tering of galaxies alone can not put stringent con- 
straints on cosmological parameters. 

These results clearly demonstrate that the observed 
SMFs at different redshifts, and the CSMFs and 2PCFs 
at low-z can all be accommodated self-consistently within 
the current ACDM model, which specifies the mass and 
bias functions of dark matter host halos, the mass func- 
tion and accretion times of dark matter subhalos, and 
the orbital decay times of subhalos within their hosts. 

6. IMPLICATIONS FOR THE ASSEMBLY OF 
GALAXIES 



The main ingredient of our model is a description of 
the relation between halo mass and the stellar mass of its 
central galaxy, as a function of redshift. In the previous 
section we have seen how this relation can be constrained 
using observational data. We now investigate what our 
models imply for the assembly of galaxies. 

6.1. The Galaxy Mass - Halo Mass Relation 

The upper panels of Fig. [131 show the M*^c - M re- 
lations at different redshifts predicted by model FIT- 
CSMF (resuhs for FIT-2PCF are fairly similar). The 
lower panels show the same relations, but this time in 
terms of the mass ratio M^,^c/M. Results are shown sep- 
arately for models that are constrained by SMFl (left- 
hand panels) and SMF2 (right-hand panels). For com- 
parison, in Fig. [141 we also plot the results obtained from 
SMFl (red dot-dashed lines) and SMF2 (red solid fines) 
at each redshift in the same panel, together with a no- 
evolution model (cyan dashed line) and two models from 
the literature (Moster et al. 2010, blue dotted lines, and 
Leauthaud et al. 2012, green long-dashed lines). 

At all redshifts the central stellar mass increases with 
halo mass. For massive halos with M > 10^'^ H'^Mq, 
the (average) central stellar mass is actually more mas- 
sive at higher redshift. Note, though, that this does not 
mean that the stellar mass of a particular central galaxy 
decreases with time. After all, the main branch mass of 
a halo also increases with time. However, it does mean 
that the central stellar mass cannot grow proportionally 
to halo mass. The straight dashed lines in Figs. [T31 and 
[131 mark a stellar mass of 10^^-^h~^ M©, which roughly 
corresponds to the largest stellar mass in the SDSS sur- 
vey (see Fig. [I9\i . If galaxies more massive than this 
are absent from the local universe, it is reasonable to 
assume that there are also no such galaxies at higher 
redshifts. Note that this is consistent with the SMFs 
shown in Figs [51 and [3 Unfortunately, since the high- 
mass ends of the high-z SMFs are uncertain (as massive 
galaxies are rare and current survey volumes are rela- 
tively small) and can be quite significantly impacted by 
the actual scatter in the stellar mass estimations (see the 
cyan, dot-dashed curves in Figs. M and [3 for an illustra- 
tion), the M*,c - M relation is not well constrained at 
the very massive end. Nevertheless, despite the uncer- 
tainties, the observational data clearly indicate that the 
central galaxy of a high-mass halo at high-z (usually too 
rare to be observed) is more massive than that of a low-z 
halo of the same mass. As we will see in more detail in 
Section 16.31 this also means that massive galaxies, cur- 
rently residing in massive halos, must have formed early 
in progenitor halos with much lower masses. The stel- 
lar masses of these massive galaxies have increased little 
during the later growth of their host halos, indicating 
that star formation in these galaxies must have ceased 
at early times, and that they have not cannibalized too 
many (massive) satellites. 

For low-mass halos, especially those with M < 
10^^ H^^Mq, there is a huge increase in the stellar mass 
for a given halo mass from z > 2 to z ~ 0. This in- 
crease is more than a factor of 10 for SMFl and a factor 
of more than 3 for SMF2 (see Fig. [TH). This has im- 
portant implications for the star formation efficiencies in 
low mass halos at high z. Afterall, the M^^^c/M ratio 
of low mass halos at low z is already far below the uni- 
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log [M, / h-'M^] log [M, / h-'M,,] 

Fig. 13.— Similar to Fig. [12] but here we show model predictions at different redshifts. Results are shown separately for SMFl (left- 
hand panels) and SMF2 (right-hand panels). For simplicity, results are only shown for model FIT-CSMF, although the results for model 
FIT-2PCF are very similar. The black dashed lines correspond to M, = 10^^'*/i~^ Mq, which roughly reflects the observed stellar mass 
limit of galaxies. 



versal baryon fraction (see lower panels of Fig. [T^ . If 
this ratio is a factor 3 to 10 lower at higher redshifts, 
it indicates that star formation in these halos must be 
extremely inefficient. As to be discussed in more de- 
tail in Section [7l this poses a serious challenge to our 
current theory of galaxy formation. Although it is stan- 
dard practice to invoke supernova feedback as a feedback 
mechanism to suppress star formation in low mass haloes, 
this process is not expected to be very efficient in haloes 
with M > 10" h-^MQ (e.g., Dekel & Silk 1986), and it 
appears that additional mechanisms are required to ex- 
plain the low inferred star formation efficiencies at high- 
z. Alternatively, one might question whether the data for 
faint galaxies is sufficiently reliable. After all, the SMFs 
at the low-mass ends obtained by D05 are significantly 
higher than those by PG08. However, as one can infer 
from a comparison between the data points and the cyan, 
dot-dashed curves in Figs. [6] and [3 the 'no-evolution' as- 
sumption (assuming the AI^^c - M relations at high-z to 



be the same as that at z = 0.1) over-predicts the high-z 
stellar mass functions at the low-mass ends by very large 
amounts, about a factor of 5 for SMFl and about a fac- 
tor of 2 for SMF2 at log(Af,/ H-^Mq) ~ 10 at redshift 
z ^ 2.0. Such large discrepancies are difficult to be ex- 
plained by observational uncertainties. Indeed, as shown 
in Marchesini et al. (2009), the SMFs of D05 at the low- 
mass ends are already the highest among a handful of 
independent observations. Hence, the implied low star 
formation efficiencies in low-mass halos at high-z have to 
be taken serious. 

6.2. Comparison with Previous Studies 

For completeness, we compare our results regarding 
the evolution of the M^,^c-M relation to a few previous 
studies. Note that in order to carry out proper compar- 
isons, one needs to take into account various systematics 
among different analyses. First, most of the analyses, ex- 
cept perhaps those based on gravitational lensing, have 
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Fig. 14. — The Mt,c - M relations in different redshift bins. In each panel we compare the predictions of model FIT-CSMF using SMFl 
and using SMF2 (shaded areas for the 68% confidence levels as indicated), respectively. For comparison, we also show in each panel the 
prediction of the "no-evolution model", which simply is the M*_c - M relation at z = 0.1 (using SMFl; dashed line), and the model 
predictions of Mostor et al. (2010; dotted line). In addition, for three redshift bins we also indicate the model predictions obtained by 
Leauthaud et al. (2012; long-dashed lines) and Wake et al. (2011; filled dots), respectively. 



to assume a cosmology, so that the result is cosmology 
dependent (see van den Bosch et al. 2003b). Second, dif- 
ferent SED (Spectral Energy Distribution) fitting codes 
may yield significantly different stellar masses even with 
the same IMF (see e.g. Behroozi et al. 2010). Finally, 
the uncertainty in the stellar masses may be redshift de- 
pendent (e.g. Drory et al. 2005). 

At low redshift, the M^,^c-M relation has been studied 
not only in the framework of HOD/CLF, but also with 
the use of weak gravitational lensing, X-ray observations, 
and satellite kinematics. As shown in More et al. (2009) 
and Leauthaud et al. (2012), these measurements are in 
general agreement with each other. As we have shown 
above, our results at redshift z = 0.1 are in good agree- 
ment with those of Yang et al. (2008), and thus are 
also in agreement with other results (see e.g.. More et al. 
2009). Here we focus on comparison of results at high 
redshift. 

Moster et al. (2010) used a method that is very similar 
to subhalo abundance matching, and stellar mass func- 
tions of Fontana et al. (2006) covering a wide range in 
redshifts, to constrain the galaxy-dark matter connection 
from z ~ 4 to the present. The resulting c - M rela- 
tions are shown as blue, dotted lines in Fig. IJH At low 
redshift, Moster et al. used the stellar mass function of 
Panter et al. (2007), which is actually very similar to the 
SMF (for SDSS DR7) we adopt (see Fig. 2). The signifi- 



cantly higher M^^^c at the massive end obtained by Moster 
et al. (2010) is likely due to the following two factors. 
First, the model of Moster et al. (2010) assumes zero 
scatter in the M*_c - M relation. As a demonstration, we 
show their model predition assuming Uc = 0.15 instead of 
fJc = in the top-left panel of Fig. [14] as the dot-dashed 
line. As one can see, assuming a finite dispersion indeed 
significantly suppresses the slope of the Af,.c - M relation 
at the high mass end. Second, Moster et al.'s results are 
based on abundance matching using a particular iV-body 
simulation of a box size of L = 100 Mpc, and may be af- 
fected by box-size effect, particularly for massive halos. 
It might be that massive halos are under-represented in 
their simulation, so that their model has to assign mas- 
sive galaxies to halos with relatively low masses. At high 
redshift, Moster et al. used a set of SMFs obtained by 
Fontana et al. (2006), where uncertainties in the esti- 
mated stellar masses might be significant. Nevertheless, 
if we use the difference in the predictions of SMFl and 
SMF2 as an indication of the observational uncertainty, 
the Moster et al. results are consistent with ours. 

We also compare our results to those of Leauthaud et 
al. (2012), at least for the three redshift bins 0.2 < z < 
0.4, 0.4 < z < 0.6, and 0.8 < z < 1.0 (green long-dashed 
curves in Fig.[T4|). Rather than using subhalo abundance 
matching, Leauthaud et al. used galaxy number densi- 
ties (SMFs), galaxy clustering and galaxy-galaxy lens- 
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ing in the COSMOS survey (Scoville et al. 2007), in or- 
der to constrain the galaxy dark matter connection from 
z ~ 0.2 to z ~ 1.0. Since they parameterized the occupa- 
tion statistics of centrals and satellites separately, they 
effectively allow for a dependence on subhalo accretion 
times, and for evolution in the satellite population after 
accretion. A comparison of their M^^c - M relations to 
ours shows that their stellar masses for central galaxies 
at a given halo mass are generally very similar to ours, 
especially for the two low redshift bins. For the highest 
redshift bin, there is noticeable discrepancy in the inter- 
mediate halo mass range (M ^ 10^^ °). Unfortunately 
the origin of the discrepancy is not clear. 

Finally, we also compare our results to those obtained 
by Wake et al. (2011) in three redshift bins: 1.0 < z < 
1.2, 1.2 < z < 1.6, and 1.6 < z < 2.0 (filled dots in Fig. 
12). These data are obtained by fitting HOD models 
to the abundance and clustering of galaxies, as function 
of stellar mass, obtained from the NEWFIRM Medium 
Band Survey (van Dokkum et al. 2009). A comparison 
with our results shows that their M^^c - M relations are 
in quite good agreement with our results. Note however, 
as pointed out in a recent paper by Leauthaud et al. 
(2011), that using different functional form for the HOD 
model might impact the average halo masses that host 
given stellar mass galaxies as well. 

6.3. The Growth of Stellar Components in Dark 
Matter Halos 

Using the redshift-dependent CSMFs obtained above, 
we can predict the growth of the stellar masses of both 
central and satellite galaxies along the main branch of 
their dark matter halos. The median stellar mass at red- 
shift z of a central galaxy that at redshift zq < z \s 
located in a host halo of mass Mq can be written as 

Af,,,(z|Mo, ^o) = M^^rXMa. z) ■ (48) 

Here, as before. Ma is the median mass at redshift z > zq 
of the main progenitor of a host halo of mass Mq at 
redshift zq. The median, total stellar mass of surviving 
satellite galaxies in the main branch can be obtained by 
integrating the CSMF of satellite galaxies: 

m,^s(z|Mo, zo) = y dlogm, TO* $s(m,|AfQ, z) . (49) 

Thus, once we know the assembly history of a dark mat- 
ter halo, it is straightforward to use our model to ob- 
tain the corresponding assembly histories of the stellar 
components of its central and satellite galaxies. In ad- 
dition, we also examine what the models predict for the 
total mass of satellite galaxies that have been accreted 
onto the main branch, which includes the surviving satel- 
lites as well as those that have been cannibalized by the 
central or disrupted by the tidal field. Note that here 
we have taken into account the evolution of the satel- 
lite galaxies after their accretion, e.g. the addtional star 
formation and passive evolution, using the parameter c. 
This component is thus given by 

m^,e.cc{z\Mo,zo) ^ [ , . /"dlogTO^m, (50) 

$e(m*|ma, Za, z) rigub (Wa , Za | Mq , Zq) . 



Here we adopt two different methods to obtain the as- 
sembly histories of dark matter halos. The first is based 
on halo merger trees extracted directly from a large A^- 
body simulation. The second uses the analytical model of 
Yll (see also Section [2^ . which relies on the fitting for- 
mula of Zhao et al. (2009) to describe the mean assembly 
history for halos of a given final mass and adopts a log- 
normal dispersion around the mean main-branch mass 
at each time. As shown in Yll, this analytical model 
accurately matches the simulation results averaged over 
halos of a given final mass. However, it does not fully 
take into account the variance among different halos, as 
it ignores the covariance of the masses at different times 
along individual main branches. Using the merger trees 
extracted from the numerical simulations does not suffer 
from this shortcoming and therefore allows us to get an 
estimate for the halo-to-halo scatter, taking full account 
of the variance in their assembly histories. 

The simulation used here is the 300/i~^Mpc simulation 
box described in Yll, which evolves 1024^^ dark matter 
particles in a WMAP5 cosmology from redshift z = 100 
to the present epoch (z = 0) using the massively paral- 
lel GADGET2 code (Springel et al. 2001; 2005). The 
particle mass and softening length are 1.87 x 10^ /i~^Mq 
and 6.75 /i~^kpc, respectively. A total of 100 outputs 
are available between z = 50 and z = 0, spaced in equal 
log(l + z) intervals. Dark matter halos were identified at 
each output using the standard friends-of-friends (FOF) 
algorithm (Davis et al. 1985) with a linking length of 0.2 
times the mean interparticle separation. Only halos with 
at least 20 particles are used. Based on halos at different 
outputs, halo merger trees were constructed. A halo in 
an earlier output is considered to be a progenitor of the 
present halo if more than half of its particles are found in 
the present halo. The main branch of a merger tree is de- 
fined to consist of all the progenitors as one climbs from 
the bottom to the top, choosing always the most massive 
branch at every branching point. These progenitors are 
the main-branch progenitors, and the time dependence of 
the main branch mass is the assembly history of the halo. 
We extracted 200 merger trees for present-day halos in 
each of the mass bins centered at log(M/ /i~^M0) = 12.0, 
13.0, 14.0 and 14.5. The median and the 68% range of 
the distribution of assembly histories thus obtained are 
indicated as solid lines and hatched areas, respectively, 
in the panels in the top row of Fig. 1151 For comparison, 
the dotted lines indicate the mean assembly histories ob- 
tained using the analytical model of Zhao et al. (2009), 
which are in excellent agreement with the solid lines. In 
general, more massive halos assemble later than smaller 
ones, which reflects the hierarchical nature of structure 
formation in the CDM model (e.g., van den Bosch 2002; 
Wechsler et al. 2002). Note that, as we have tested our 
results using another set of halo merger trees constructed 
with a different method from simulations of the WMAPl 
cosmology (Han et al. 2011), we found that our results 
are robust against the change in cosmology and in the 
method of merger-tree construction. 

In the following subsections, we use these halo assem- 
bly histories to examine in detail how the different stellar 
components in dark matter halos (central galaxies, satel- 
lite galaxies, and halo stars) grow with time. 

6.3.1. Central Galaxies 
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Fig. 15. — The growth of host (main branch) halo mass (top row panels), the stellar mass of central galaxies (second row), the total stellar 
mass in surviving satellite galaxies (third row), the total stellar mass brought into the main branch by sub-halos (fourth row), and the 
stellar mass of the "disrupted" component (bottom row, see text for definition), all as function of redshift. These results are obtained using 
the best-fit parameters of model ID=10, listed in Table |3] (using WMAP5 cosmology and SMFl). Results in different columns correspond 
to host halos with masses centered at log(M//i~^M0) = 14.5, 14.0, 13.0 and 12.0, from left to right. The solid line and shaded area in each 
panel indicate the median and 68% percentile obtained from 200 simulated halo merger trees. The dotted line is the prediction obtained 
using the analytical model of Yll. The two horizontal dotted lines in the second row indicate the maximum and half-maximum values of 
the median stellar mass of the central galaxies. The dashed lines in the third, fourth and bottom rows are the same as the solid lines in 
the top panels. Finally, for comparison, the dot-dashed lines (in all rows except the top one) are model predictions obtained using the best 
fit parameters of model ID=10 listed in Table |4] (using WMAP5 cosmology and SMF2). 



Using model FIT-CSMF for the WMAP5 cosmology 
(ID=10 in Tables [3] and HI , we populate the host halos in 
the simulation outputs with central galaxiefl Here we 
ignore the scatter in the central stellar mass - halo mass 
relation, simply because it is relatively small and because 
it is difficult to take it into account properly owing to pos- 
sible covariances between the stellar masses in the main 

5 We have also used the FIT-2PCF models (ID=9) instead, and 
found very similar results. For the sake of brevity, we will only 
show results for model FIT-CSMF. 



branch of a halo at different times. The results are shown 
in the second row of Fig. I15[ where different panels corre- 
spond to different present-day halo masses, as indicated 
in the top panels. The solid lines and shaded areas in- 
dicate the median and 68 percentiles of M^,,c{z\Mo, zq) 
obtained from 200 halo merger trees. For comparison, 
the red, dotted curves show the M^,^c{z\Mo, zq) obtained 
from Eq. (|48| . Note that these agree extremely well with 
the results obtained from the merger trees in the numer- 
ical simulations. 
In each panel, the two horizontal doted lines indicate 
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1+z l+z 1+z 1+z 

Fig. 16. — Similar to the lower four rows of Fig. 1151 but this time we h ave normalized all stellar masses by the median, instantaneous 
main-branch halo mass, (the solid lines in the top panels of Fig. I15II : i.e., the various lines indicate the median stellar-to-halo mass 
ratios (in units of h~^). From the top row to the bottom row, the results are shown for the stellar mass in central galaxies (first row), the 
total stellar mass in surviving satellite galaxies (second row), the total stellar mass brought into the main branch by sub-halos (third row), 
and the stellar mass of the "disrupted" component (bottom row). Line styles are the same as in Fig. 1151 



the maximum and half of the maximum of the median 
stellar mass of the central galaxies at z = 0. These can 
be used to read off when the central galaxies, on average, 
acquired half their present day stellar masses; for centrals 
in present-day halos of log(M//i-iMQ) - 12.0, 13.0, 14.0 
and 14.5 this occurs at z ^ 0.5, 1.0, 1.5 and 2.0, respec- 
tively. Thus, although more massive halos assembe their 
masses later, their central galaxies actually assembled 
earlier than centrals in lower-mass halos (see also Con- 
roy & Wechsler 2009). Note how the stellar mass assem- 
bly histories of central galaxies reach a plateau in halos 
with M > 10^^ h~^MQ, indicating that central galaxies 
in massive halos grow very little in stellar mass at low 
redshifts. In particular, central galaxies in clusters with 
M ^ lO^'*'^/i"^M0 have not grown in mass (significantly) 
since z ~ 1.0. In contrast, central galaxies in Milky- Way 
sized halos with M ~ 10^^ h~-^MQ continue to grow in 
stellar mass down to the present day. The slow growth of 
central galaxies in massive halos suggests not only that 
their star formation must have been quenched by some 
processes, but also that their increase in stellar mass due 
to the accretion of satellite galaxies cannot have been 
signficant. As discussed in Section [3 this has important 



implications for our understanding of galaxy formation 
and evolution. 

Another interesting way to depict the growth of the 
stellar mass of central galaxies is to show the ratio be- 
tween stellar mass and host halo mass (the main branch 
mass) at the redshift in question. These curves are 
shown in the top panels of Fig. [16] using the same line 
styles as in Fig. [151 As one can see, central galaxies 
in host halos of different final masses reach a maximum 
of about (0.015 — 0.03)/i~^, quite independent of the fi- 
nal host halo mass. However, the maximum is reached 
at very different redshifts of z ~ 5, 4, 1.5 and for 
log(Af//i-iM0) - 14.5, 14.0, 13.0, and 12.0, respectively. 
Interestingly, all these redshifts roughly correspond to a 
similar main branch mass of ~ 10^^ h~^MQ. 

Finally, in order to illustrate the effects of the uncer- 
tainties in the high-z SMF measurements, the magenta 
dot-dashed curves in Figs. [T5| and \W\ show the median 
stellar masses obtained using the D05 SMFs (case SMF2, 
the best-fit parameters of which are listed under ID=10 
in Table [U . The results are qualitatively the same as 
those obtained from SMFl. However, the stellar masses 
of central galaxies in low-mass halos at high-z are higher 
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than those obtained from SMFl. This simply reflects 
that the high-z SMFs of DOS are significantly steeper at 
the low mass end than the PG08 SMFs (see Fig.|l|). 

6.3.2. Satellite Galaxies 

Next we use the merger trees in our simulation box 
to investigate how the population of surviving satellite 
galaxies grows in stellar mass. The results are shown 
in the third row of panels in Fig. [15] Once again the 
solid lines and shaded areas indicate the median and 
68 percentiles of m^,^s{z\Mo, zq) obtained from 200 halo 
merger trees, while the red dotted lines indicate the 
(m*,s)(z|Mo, zo) obtained from Eq. (|^ using the an- 
alytical model of Yll for the halo assembly history. The 
black, dashed lines indicate the median of the stellar mass 
of central galaxies, taken from the panels in the second 
row, and are shown for comparison. As one can see, the 
total mass of satellites only becomes larger than that of 
the central galaxy in massive halos at low redshifts. In 
Milky- Way sized halos the total satellite mass never ex- 
ceeds ~ 10% of the central mass. For all halos, satellite 
galaxies contain less than a few percent of the total stellar 
mass at all z > 3. The total satellite mass in a massive 
halo continues to grow all the way to the present time, 
whereas the central galaxy ceases to grow once the halo 
mass reaches ~ 10^^ /i~^Mq. As shown by the magenta, 
dot-dashed lines, the use of SMF2 predicts higher satel- 
lite masses at high-z than the use of SMFl, especially 
in low-mass halos. Once again this reflects the fact that 
the high-z SMFs of D05 are steeper at the low mass end 
than their PG08 counterparts. 

The ratios between the total satellite mass and the 
main branch halo mass are shown in the second row of 
Fig. mi For low- mass halos and for massive halos at early 
times, this ratio increases rapidly with time, implying 
that the total satellite mass grows faster than the halo 
mass. The reason for this is that for a given halo mass, 
the stellar mass of central galaxies in low-mass halos in- 
creases rapidly with time (see upper panels of Fig. [16]), 
and so the ratio between the stellar mass and dark mat- 
ter mass of the accreted matter increases rapidly with 
time. For massive halos, the ratio flattens at low-z to a 
value ^ 0.002/i^^, corresponding to a halo mass to stel- 
lar mass ratio of ~ 200ft. which is close to the universal 
mass-to-light ratio. 

6.3.3. Halo Stars and Star Formation in Central Galaxies 

The third component we consider is the total stellar 
mass that is brought into the main branch by all satel- 
lite galaxies at the time of their accretion. Here we 
have taken into account the later evolution effect of the 
satellite galaxies, such as the addtional star formation 
and passive evolution. This mass contains not only the 
mass of the 'surviving' satellite galaxies, but also that of 
the satellites that have been cannibalized or disrupted. 
The mass growth histories of this mass component are 
shown in the fourth row of Fig. I15[ where as before 
the solid lines and shaded areas depict the medians and 
68 percentiles obtained from 200 merger trees, respec- 
tively. The dotted lines show the results obtained by 
integrating Eq. (jSO)) . The corresponding ratios between 
this mass component and the main branch halo mass are 
shown in the third row of panels of Fig. 1161 In halos 
with M < 10" /i^^Mq, the total accreted steUar mass 



is roughly similar to that of the surviving satellite pop- 
ulation. In halos with M > 10^^ h~^MQ, however, the 
total accreted stellar masses are significantly larger than 
those in the surviving satellite population, especially at 
low rcdshift (z < 1.0). This is more clearly seen in the 
bottom rows of Fig. [15] and [TH] which show the differ- 
ences between the median accreted mass and the median 
mass of the surviving satellites. This indicates that in 
massive halos a significant fraction of the total stellar 
mass brought in by subhalos is actually not associated 
with surviving satellite galaxies and must exist in other 
forms. For convenience, we will refer to this stellar mass 
component as the "disrupted" component. 

There are two possible fates for the disrupted compo- 
nent: either it is accreted by the central galaxy or it gives 
rise to a population of halo stars (also referred to as in- 
tracluster stars, or ICS, in the case of massive halos). 
Unfortunately, our model cannot uniquely discriminate 
between these two components, as the mass of the central 
galaxy can also grow through star formation. However, 
we can obtain important constraints by comparing the 
growth of the disrupted component to that of the cen- 
tral galaxies (black, dashed lines). This shows that cen- 
tral galaxies in halos with M < 10^'^ h~^MQ must have 
grown mainly through in situ star formation, simply be- 
cause the total mass in the disrupted component is much 
smaller than that of the central. For more massive halos, 
however, the situation could be very different. Here the 
total stellar mass in the disrupted component is actually 
larger than that of the central galaxy. Hence, the exis- 
tence of ICS is inevitable; the total stellar mass locked up 
in the stellar halo can be more than 2 times as massive 
as the central galaxy at redshifts z < 1.5 (see also Pur- 
cell et al. 2007). As one can see from the left two panels 
in the bottom row of Fig. [T5] although the total mass 
of the disrupted fraction is sufficient to form the cen- 
tral galaxy in massive halos, in reality central galaxies in 
these halos must have acquired their stellar mass mainly 
through in situ star formation, rather than through the 
accretion of satellites (see also results along this line ob- 
tained from hydro-dynamical simulations by e.g., Naab 
et al. 2009; Oser et al. 2010; 2012; Hrischmann et al. 
2012). This can be inferred from the fact that these cen- 
tral galaxies have already build up most of their stellar 
mass at high redshifts, when the mass in the disrupted 
component is too small to make a dominant contribution. 
Accretion can only be significant at z < 2 in massive ha- 
los, when the disrupted component is sufficiently massive 
compared to the mass of the central. However, at these 
low redshifts, the central galaxies in these massive halos 
no longer seem to grow much in stellar mass. Hence, 
even central galaxies in massive halos must have mainly 
grown via in situ star formation. 

Finally, we point out that the use of SMF2 leads to 
masses in the disrupted component that are quite simi- 
lar to those obtained using SMFl. The higher total ac- 
creted mass predicted with SMF2 for low-mass halos at 
high-z is largely due to the higher mass in the survived 
satellites. Note also that the analytical model predicts 
larger disrupted components in low-mass halos at high-z 
than using the merger trees in the numerical simulation 
(cf solid black lines and dotted red lines in the bottom 
panels of Figs. [15] and [TO]) . This is mainly due to the 
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Fig. 17. — The projected 2PCFs for galaxies in different stellar mass and redshift bins. Here we compare model predictions (various 
lines) with the observational measurements obtained by Meneux et al. (2009) from zCOSMOS (open circles with error- bars). The long- 
dashed, dotted, and dashed lines show the predictions of models FITO, FITl, and FIT-CSMF, respectively, all of which have used SMFl 
as constraints for the stellar mass functions. For comparison, the solid lines show the predictions for model FIT-CSMF using SMF2. 
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Fig. 18. — Predictions for the projected 2PCFs in different ste llar mass bins at redshifts z = 2.0 (upper panels) and z = 3.0 (lower 
panels). Results are shown for the same four models as in Fig. 1171 using the same line-styles, as indicated. 
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fact that small progenitors are not well resolved in the 
simulation at high-z (see the top-right panel of Fig. [15]). 

6.4. Galaxy Clustering at High Redshift 

With the model described in Section 13. 2[ we are also 
able to predict the correlation functions of galaxies at 
high-z. In Fig. [T7] we show the model predictions for 
FITO, FITl, and FIT-CSMF using SMFl as well as for 
FIT-CSMF using SMF2 (the predictions of FIT-2PCF 
are similar to those of FIT-CSMF). The model pre- 
dictions are compared with the observational results of 
Meneux et al. (2009), which are shown as open circles 
with errorbars. As one can see, all our models under- 
predict the correlation function on large scales, particu- 
larly for galaxies in the redshift range 0.5 < z < 0.8. As 
pointed out in Meneux et al. (2009), and as briefly men- 
tioned in Sect ion l6 . 1 1 above . these enhanced correlations 
on large scales are likely due to the fact that the zCOS- 
MOS field happens to correspond to a high density region 
of the Universe. Nevertheless, the results show that dif- 
ferent models predict very different correlations on small 
scale (rp < l/i~^Mpc), once again demonstrating that ac- 
curate clustering measurements on small scales can pro- 
vide important constraints on the galaxy-dark matter 
connection (in particular regarding the dynamical evo- 
lution of satellite galaxies) . 

In Fig.[T8]we show model predictions for the projected 
2PCFs at z = 2.0 (upper panels) and z = 3.0 (lower 
panels). Here again results are shown for FITO, FITl, 
and FIT-CSMF using SMFl, as weh as FIT-CSMF us- 
ing SMF2. Different panels correspond to galaxies in dif- 
ferent stellar mass ranges, as indicated. The projected 
2PCFs on large scales decrease with both increasing red- 
shift and decreasing stellar mass. Note that the differ- 
ences among the different models are extremely small at 
large scale (r-p > 1 ft,~^Mpc). On small scales our models 
predict that the 2PCFs are greatly enhanced relative to 
a simple extrapolation from large scales, except for the 
FITO models which has zero satellites, and therefore no 
one-halo term. Note that this enhancement is stronger 
for more massive galaxies, and that it is fairly different 
for the different models. Hence, as for the low-z 2PCFs, 
the clustering strength on large scales does not provide 
much more information than the SMFs themselves (see 
also Moster et al. 2010), while that on small scales does. 

7. SUMMARY AND DISCUSSION 

In this paper we have developed a new and self- 
consistent model that describes the galaxy-dark matter 
connection over cosmic history. Unlike the popular abun- 
dance matching technique, our model takes account of 
the fact that sub-halos in a host halo are accreted at dif- 
ferent times, so that the stellar masses associated with 
them are likely to depend on both their halo masses at ac- 
cretion and their accretion times. In addition, our model 
allows for dynamical evolution of the satellite population 
(mass stripping, tidal disruption, cannibalism). 

We have used galaxy stellar mass functions observed 
in the redshift range < z < 4, together with the con- 
ditional stellar mass function and two-point correlation 
functions (both at z ~ 0.1, and both obtained from 
SDSS), to constrain the evolution of the galaxy-dark 
matter connection from z ~ 4 to the present. The re- 
lation between halo mass and stellar mass thus obtained 



is used, together with simulated and theoretical halo as- 
sembly histories, to predict how the masses of different 
stellar components in dark matter halos (central galax- 
ies, satellite galaxies, and halo stars) evolve as function 
of time. We also used our model to predict the 2PCFs of 
high-z galaxies as a function of their stellar masses. Our 
main findings can be summarized as follows. 

1. Our model provides a reasonable fit to all the data 
within the observational uncertainties, which indi- 
cates that the current ACDM model is consistent 
with a wide variety of data regarding the galaxy 
population across cosmic time. 

2. At low-z, the stellar mass of central galaxies, M*^c, 
increases with halo mass as Af" '^ and Af~^ ° at the 
massive and low-mass ends, respectively. The ra- 
tio M^^^jM reveals a maximum of ~ 0.03 at a halo 
mass M ~ 10"-* /i"^Mq. The fact that this maxi- 
mum value is much lower than the universal baryon 
fraction (^ 0.17) reflects the overall inefficiency of 
star formation. At higher redshifts the maximum 
in Mif^c/M remains close to ~ 0.03, but shifts to 
higher halo mass (see also Wake et al. 2011). 

3. The Af*^c-Af relations obtained here are roughly 
in agreement with previous studies (e.g., Moster 
et al. 2010; Wang & Jing 2010; Leauthaud et 
al. 2012). However, our model allows us to put 
constraints on the satellite population so that we 
can interpret the results in terms of the evolution of 
galaxies in dark halos. We find that low-mass satel- 
lite galaxies can significantly increase their stellar 
masses after accretions by their host halos. 

4. The time-scale for the disruption of satellite galax- 
ies is about the same as the dynamical friction time 
scale of their subhalos obtained in A^-body sim- 
ulations, suggesting that most disruption occurs 
once the galaxy has 'depleted' its orbital energy, 
at which point it is either accreted by the central 
galaxy or torn apart by the strong tidal field. 

5. The stellar mass assembly history for central galax- 
ies is completely decoupled from the assembly his- 
tory of its host halo. Initially the ratio M*^c/-W 
increases rapidly with time, until the halo mass 
reaches ~ 10^^ /i~^M0, at which point M^^c/M ^ 
0.03. Once M > lO^^ft-^M©, there is little growth 
in M^,^c, causing the ratio M^,,c/M to decline. 

6. Most massive centrals assemble their stellar mass 
earlier than their less massive counterparts, even 
though their host haloes assemble later. In Milky- 
Way sized halos more than half of the central stellar 
mass is assembled at z < 0.5, while brightest clus- 
ter galaxies had already grown to half their present 
day stellar mass by z ~ 2.0. 

7. The accretion of satellite galaxies contributes little 
to the formation of central galaxies in small ha- 
los. Hence, most of their stars must have formed 
in situ so that their assembly histories closely re- 
semble their star formation histories. In massive 
halos more than a half of the stellar mass of the 
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central galaxy has to be formed in situ at z > 2. 
The accretion of stars in satellites by the centrals 
only becomes significant at z < 2, when the total 
mass available for the accretion reaches a signifi- 
cant fraction of the mass of the central galaxies. 

8. In massive halos, the total mass in halo stars is 
more than two times that of the central galaxy, 
while in Milky- Way sized halos the mass in halo 
stars is much smaller ( < 10%) than that of the 
central galaxy. 

9. The 2PCFs of galaxies on large scales do not pro- 
vide much more constraints than those already pro- 
vided by the stellar mass functions of galaxies (see 
also Moster et al. 2010). On small scales, however, 
the 2PCFs provides important constraint on how 
satellite galaxies evolve within dark matter halos. 
Future observations of the small scale clustering 
strength of galaxies at high-z will therefore be im- 
portant for our understanding of galaxy formation 
and evolution. Our models predict that the high-z 
2PCFs on small scales are much steeper than at 
low-z, especially for more massive galaxies. 

These results have important implications for our un- 
derstanding of galaxy formation and evolution. Here we 
highlight a few issues which we believe to be particularly 
interesting. 

One of the implications of our results is that star for- 
mation in low mass halos at high-z must be extremely in- 
effecient. This basically follows from the strong evolution 
in the M^,_c-M relation that we infer between high-z and 
the present at the low mass end: halos of 10^^ h~^MQ 
host central galaxies at z ^ A that are orders of mag- 
nitude less massive than centrals in halos of the same 
mass at z = (see Fig. [T3|) . Note that even at z = 
central galaxies in halos with M — 10^^ H^^Mq have 
a stellar mass that is only 2-3 percent of the universal 
baryon fraction. Hence, only a minute fraction of the 
baryons associated with low mass halos at high-z are 
turned into stars. It is typically assumed that star for- 
mation is suppressed in low-mass halos due to supernova 
feedback and/or photo-ionization. However, it has be- 
come clear in recent years that either the standard treat- 
ment of these processes is inadequate, or that additional 
mechanisms are needed to help suppress star formation 
in low mass halos. In particular, recent observations have 
indicated that the specific star formation rates (SSFRs) 
of star forming galaxies, at fixed stellar mass, increase 
with redshift from z = until they reach a plateau at 
z>2 (see e.g.. Stark et al. 2009; Labbe et al 2010a,b; 
Gonzalez et al. 2010). As pointed out by several studies, 
this constancy of the SSFRs is inconsistent with predic- 
tions of 'standard' models for galaxy formation, which 
instead predict SSFRs that continue to increase with red- 
shift (Lo Faro et al. 2009; Dutton, van den Bosch & Dekel 
2010; Bouche et al. 2010; Firmani & Avila-Reese 2010; 
Wcinmann, Neistein & Dekel 2011; Lacey et al. 2011). 
Simply increasing the supernova feedback efficiency in 
order to lower the star formation rates at high-z is not 
a solution, because it has no impact on the SSFRs, as it 
changes both the SFR and the stellar mass (Dutton et 
al. 2010; Bouche et al. 2010). In addition, it will result in 



too little stellar mass at low and intermediate redshifts, 
not to mention the fact that the supernova feedback ef- 
ficiencies that are typically invoked are already unrealis- 
tically high (e.g., Benson et al. 2003). What seems to be 
needed instead is a modification of the 'standard' model 
that (i) causes a larger suppression of star formation in 
low mass halos at high-z, and (ii) a boost in star forma- 
tion at intermediate redshifts (z ~ 1 — 2). Although there 
has been no shortage of ideas (e.g.. Mo et al. 2005; Lo 
Faro et al. 2009; Bouche et al. 2010; Lacey et al. 2011; 
Krumholz & Dekel 2011; Weinmann, Neistein & Dekel 
2011; Wang, Weinmann & Neistein 2011; Avila-Reese et 
al. 2011), it is clear that much more work is required 
before this outstanding problem in galaxy formation and 
evolution is adequately addressed. We believe that our 
constraints on the stellar assembly histories of (central) 
galaxies may play an important role in testing and cal- 
ibrating these new models. In particular, the rapid in- 
crease in M^^c/M at high-z (before the host halo reaches 
a mass M ~ 10^^ /i~^Mq) indicates that the specific star 
formation rate of (central) galaxies is much higher than 
the specific growth rate of its dark matter halo, which 
may hold important clues. 

Another intriguing feature from the stellar mass assem- 
bly histories is the fact that M,_c/-M reaches a maximum 
of 0.03 once M reaches a mass of ~ 10^^ h'^MQ. This 
feature seems to hold independent of the final host halo 
mass of the central galaxy. Since, as we have shown, the 
stellar mass growth in host halos with M < 10^^ h~^MQ 
is dominated by in situ star formation (rather than ac- 
cretion), this suggests that something quenches star for- 
mation in central galaxies once its halo mass reaches 
lO^^ft-~^M0. Interestingly, this mass scale is very similar 
to the cold-mode to hot-mode transition scale (Birnboim 
& Dekel 2003; Keres et al. 2005), suggesting that the 
quenching of central galaxies coincides with the forma- 
tion of a hot gaseous halo. This is indeed what seems 
to be required in order to explain the observed bimodal- 
ity in colors and/or specific star formation rates of the 
galaxy population (e.g., Cattaneo et al. 2006; Birnboim 
et al. 2007). Different mechanisms have been invoked 
to explain why this transition is associated with a shut- 
down of star formation in the central galaxies, ranging 
from AGN feedback (e.g., Tabor & Binney 1993; Ciotti & 
Ostriker 1997; Croton et al. 2006; Bower et al. 2006; Hop- 
kins et al. 2006) and gravitational heating (e.g., Fabian 
2003; Khochfar & Ostriker 2008; Dekel & Birnboim 2008; 
Birnboim & Dekel 2011) to thermal conduction (e.g., 
Kim & Narayan 2003) and turbulence (e.g., Zhu, Feng 
& Fang 2011). Although our results are unable to pro- 
vide direct insight into the exact quenching mechanism, 
it is fascinating (or at least reassuring) that our analysis, 
which uses no (direct) data on star formation rates or the 
color-bimodality of the galaxy population, also comes to 
the conclusion that galaxies are quenched once their halo 
mass reaches ^10^^ h~^MQ. 

Another important result regards the implied build- 
up of stellar halos. We have argued in Section 16.3.31 that 
central galaxies build up the vast majority of their stellar 
mass via in situ star formation (the only possible excep- 
tion are the most massive centrals, which may have build 
up as much as roughly half their stellar mass via accre- 
tion). This implies that the difference between the total 
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mass in accreted satellites and surviving satellites, what 
we called the "disrupted" component in Section I6.3.3[ 
must have given rise to a stellar halo. This result is 
consistent with a number of recent studies which have 
argued that reconciling halo occupation statistics with 
halo merger rates requires that a significant fraction of 
satellite galaxies is indeed tidally disrupted (e.g., Conroy, 
Ho & White 2007; Conroy, Wechsler & Kravtsov 2007; 
Kang & van den Bosch 2008; Yang, Mo & van den Bosch 
2009a). As these studies have argued, satellite disruption 
is an important ingredient of galaxy formation, which has 
hitherto been largely ignored in semi-analytical models. 
Properly accounting for satellite disruption may allevi- 
ate the problem with the excessive growth of massive 
galaxies (e.g., Monaco et al. 2006; Conroy, Wechsler & 
Kravtsov 2007, Brown et al. 2008), with the overabun- 
dance of satellite galaxies in the semi- analytical models 
(Liu et al. 2010), and may even be important for under- 
standing the observed metallicities of satellite galaxies 
(Pasquali et al. 2010). Note that our model predicts that 
the fraction of stars associated with the stellar halo is an 
increasing function of halo mass (see also Purcell, Bul- 
lock & Zentner 2007 and Henriques, Bertone & Thomas 
2008), in qualitative agreement with observations (Gon- 
zalez, Zaritsky & Zabludoff 2007). 

Finally we point out that our results may have im- 
portant implications for the 'missing satellite problem', 
which refers to the fact that the predicted subhalo count 
in a Milky- Way (MW) sized halo vastly exceeds the num- 
ber of observed satellite galaxies in the MW (Moore 
et al. 1999; Klypin et al. 1999). Ahhough theoretical 
progress combined with the discovery of a large popula- 
tion of new (typically ultra-faint) satellite galaxies in the 
Milky Way (e.g., Belokurov et al. 2007) has at least par- 
tially alleviated this concern about a mismatch between 
the numbers of low-mass subhalos and faint MW satel- 
lites, it has become clear that it is difficult to reconcile 



the observed abundance of satellite galaxies with their 
kinematics: whereas the line-of-sight velocity dispersions 
of the satellites suggest that they reside in subhalos of 
relatively low mass, their abundance suggests they reside 
in much more massive subhalos (e.g., Madau, Diemand 
& Kuhlen 2008; Boylan-Kolchin, Bullock & Kaplinghat 
2011). This problem is particularly acute if one uses 
subhalo abundance matching techniques to link satellite 
galaxies to subhalos (e.g., Boylan-Kolchin et al. 2011, 
see also Busha et al. 2011). As we have argued, subhalo 
abundance matching is not self-consistent in that it does 
not account for dependence on subhalo accretion times. 
Since our model shows that there is very strong redshift 
dependence in the M^:^c-M relation at the low mass end, 
with low mass halos at higher z hosting central galaxies 
that are much less massive, the large scatter in subhalo 
accretion times implies a huge amount of scatter in the 
ratio between the stellar mass of satellite galaxies and 
their halo mass at infall. Such a large scatter may be 
exactly what is needed to reconcile the observed satel- 
lite population of the MW with predictions from ACDM 
cosmologies (e.g,., Madau et al. 2008), and may already 
have observational support from the fact that the ultra- 
faint dwarfs appear to reside in halos of similar mass 
as the much more massive 'classical' dwarfs (Strigari et 
al. 2008; Walker et al. 2009). We intend to return to this 
intriguing issue in more detail in a future paper. 
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APPENDIX 

A. THE PROJECTED TWO-POINT CORRELATION FUNCTION OF GALAXIES AND ITS 
DEPENDENCE ON STELLAR MASS FROM SDSS DR7 

In this Appendix we present our measurements of the projected 2PCFs for galaxies in different stellar mass bins. We 
use the New York University Value-Added Galaxy Catalogue (NYU-VAGG; Blanton et al. 2005), which is based on 
SDSS DR7 (Abazajian et al. 2009) but with an independent set of significantly improved reductions. From the NYU- 
VAGC, we select all galaxies in the Main Galaxy Sample with an extinction corrected apparent magnitude brighter 
than r — 17.72, with redshifts in the range 0.01 < z < 0.20 and with a redshift completeness Cz > 0.7. This gives a 
sample of 639, 359 galaxies with a sky coverage of 7748 square degrees. For each galaxy, we estimate its stellar mass 
using the fitting formula of Bell et al. (2003). 

A.l. Galaxy samples 

Galaxies are separated into five stellar mass bins, with \og[Mt/ h'^MQ] = [9.0, 9.5], [9.5, 10.0], [10.0, 10.5], [10.5, 11.0] 
and [11.0, 11.5], respectively. These samples are referred to as 'Mass-limited' samples, and are indicated in the stellar 
mass versus redshift plot shown in the upper left-hand panel of Fig. 1191 The number densities of galaxies in the 
five stellar mass bins as function of redshift are shown as histograms in the upper right-hand panel of Fig. [T21 For 
clarity the number densities are scaled by constant factors. The advantage of using these 'Mass-limited' samples 
is that it maximizes the number of galaxies in each mass bin. However, the drawback is that these samples are 
not homogeneous in the radial direction, and therefore not straightforward to compare with model predictions. We 
therefore also construct another set of samples as follows: For each of the five stellar mass bins, we first find the 
maximum number density of galaxies (at some redshift). Next we determine the minimum and maximum redshifts at 
which the number density of galaxies (in that stellar mass bin) drops to half of the maximum value. Only galaxies 
within these two redshift limits are kept in the sample. In what follows we refer to the five samples thus selected as 
'Volumel' samples, and the corresponding number densities as function of redshift are shown as the shaded areas in 
the upper right-hand panel of Fig. ]19\ 



The Galaxy - Dark Matter Connection 



31 



12 



11 



10 



O 



1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 










- .Ill 


I'I'-ll'J 

1 


19 
l<l 

( 


m 

1 








i>,->-..:.N><>.l.»l'i|itlltllllllilllllllli 


IlilllMi 




iiiiiiicii 


,.J 




::■■;'!:?! 

- 


' ''||P 




lil'jll-liiii- 














1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 ~ 



0.05 0.1 0.15 0.2 
z 




0.15 0.2 



12 



11 

N 



10 



00 

o 



1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 






■iiyiiiiiii'!™ 










: 

ill 










~i/iiiiiiiiiiiiiiiiii~ 



0.05 



0.1 
z 



0.15 



0.2 




0.2 



Fig. 19. — The distributions of stellar mass versus redshift for galaxies in the SDSS DR7 (left-hand panels), and the number density 
distribution of these galaxies as a function of redshift (right-hand panels). Galaxies within five stellar mass bins are selected. The selection 
criteria for these five samples are specified using their redshift distributions (histograms and/or the shaded areas). See text for details 
regarding the constructions of these samples. 



Since the SDSS galaxy sample is flux limited, red and blue galaxies of the same stellar mass suffer from different 
selection effects because of their different stellar mass-to-light ratios. Using conservative limits, van den Bosch ct 
al. (2008) have shown that the stellar mass completeness limit as a function of redshift for the SDSS catalogue is given 
by 

4.852 4- 2.24:6\og Dl{z) + 1.123 log(l + z) - 1.186z 



log[M,,uJ{h-' Mq)] 



1 - 0.067Z 



(Al) 



Using this limit we also construct five samples in stellar mass that are 'complete'. The corresponding selection criteria 
are shown in the lower left-hand panel of Fig. [T21 while the shaded areas in the lower right-hand panel indicate the 
corresponding number densities as function of redshift. In what follows we refer to these five samples as 'Volume2' 
samples. 

To estimate the two-point correlation functions, one needs to construct random samples to normalize the galaxy-pair 
counts (see below). In order to take account of the overall selection effects in the SDSS catalogue, we first construct a 
random sample using the SDSS luminosity function and then apply various observational selection effects (magnitude 
limits, redshift completeness, sky boundary) according to the SDSS survey mask. Since galaxies in different stellar 
mass bin do not follow the overall luminosity function, we re-assign a redshift to each mock galaxy by randomly 
sampling the redshifts of the SDSS galaxies in the corresponding stellar mass bin. Hence, the random samples have 
exactly the same redshift distributions as their corresponding SDSS samples. By applying the same redshift cuts, we 
obtain the random samples for the 'Mass-limited', 'Volumel' and 'Volume2' galaxy samples, respectively. 
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Fig. 20. — The projected 2PCFs of galaxies in different stellar mass bins, as indicated by the log[A/»/(/i~^ M0)]-values in brackets. 
Open circles with error bars are the results for the Volumel samples. Solid and dashed lines show the results for the Volume2 samples and 
Mass-limited samples, respectively. 



TABLE 5 

The projected 2PCFs measured from the SDSS DR7 within different stellar mass bins. 



ID 












Wpirp) 




Wp(rp) 




u.p(rp) 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 






[9.0, 9.5] 


[9.5, 10.0] 


[10.0, 10.5] 


[10.5, 11.0] 


[11.0,11.5] 


1 


0.12 


373.399 


200.617 


388.815 


83.198 


393.234 


30.913 


597.466 


90.192 


2423.422 


1684.429 


2 


0.19 


259.170 


92.989 


288.384 


55.778 


287.154 


25.668 


411.860 


24.702 


1340.480 


391.143 


3 


0.30 


221.851 


69.188 


228.047 


40.157 


223.528 


14.039 


267.002 


23.886 


746.369 


123.150 


4 


0.48 


152.048 


46.861 


161.172 


25.446 


164.613 


7.819 


177.249 


10.115 


366.864 


30.737 


5 


0.75 


112.231 


32.711 


111.791 


14.929 


113.182 


3.683 


129.185 


5.566 


241.362 


21.748 


6 


1.20 


67.752 


12.318 


69.407 


7.387 


73.073 


1.938 


83.692 


3.573 


137.289 


7.388 


7 


1.90 


42.485 


6.752 


45.751 


3.630 


49.318 


0.960 


57.954 


2.935 


98.301 


3.801 


8 


3.01 


27.063 


2.678 


30.453 


2.027 


34.459 


0.776 


41.691 


1.673 


66.697 


2.229 


9 


4.77 


16.820 


1.987 


20.490 


0.995 


24.065 


0.617 


29.157 


1.126 


45.976 


1.532 


10 


7.55 


10.066 


1.942 


12.326 


0.888 


15.995 


0.983 


20.199 


1.041 


31.164 


0.951 


11 


11.97 


5.762 


0.959 


6.884 


0.677 


9.983 


0.895 


12.782 


0.942 


18.996 


0.569 


12 


18.97 


3.637 


1.183 


4.021 


0.539 


5.171 


0.673 


7.316 


0.618 


10.485 


0.414 


13 


30.07 


2.613 


0.369 


2.850 


0.229 


3.470 


0.229 


4.050 


0.647 


5.545 


0.297 


14 


47.66 


0.486 


0.448 


0.532 


0.463 


1.306 


0.922 


2.619 


0.762 


2.473 


0.213 



Note. — Column (1): a counting ID. Column (2): the projected conioving distances in unit of h ^Mpc. Column 
(3-12): the projected 2PCFs,u;p(rp), and their errors, cr^up , for galaxies in different stellar mass bins, as indicated (in 

terms of log(M,/ /i'^Mq)). 



A. 2. The estimator of the two-point correlation function and the results 
We estimate the two-point correlation function (2PCF), £,{rp,r.^), for galaxies in each sample using the following 



estimator, 
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where (DD), (-R-R), and (DR) are, respectively, the number of galaxy-galaxy, random-random, and galaxy-random 
pairs with separation (rp,r7r) (Hamilton 1993). The variables Vp and r^r are the pair separations perpendicular and 
parallel to the line-of-sight, respectively. Explicitly, for a pair of galaxies, one located at Si and the other at S2 with 
Si — cziVi/Ho, we define 

s • 1 



111 



S • S — TT^ 



(A3) 



where 1 = (si -|- S2)/2 is the line of sight intersecting the pair, and s = Si — S2. Since redshift-space distortions only 
affect Ttt, the projection of f (r-p, r^) along the axis is not sensitive to peculiar velocities, and directly related to the 
real-space correlation function. This projected 2PCF, u;p(rp), is estimated using 



^(rp,r^)dr^ 



2E^(^P' 



r) A'V 



(A4) 



(Davis & Peebles 1983). In our analysis, the summation is made over fc = 1 to 40, which, for our adopted bin width 
of ArTT — 1 h^^Mpc, corresponds to an integration from — h^^Mpc to jv = 40 h~^Mpc. In order to obtain the 
error bars and the covariance matrix for the projected 2PCFs, we use 200 bootstrap resamplings of the galaxies in 
consideration. The covariance matrix is obtained using 



Cij = Cov(t 



1 



N 



N 



1 ^ 

1=1 



I 



(A5) 



200, and ^ represents the value of the projected 2PCF of the ith bin in the Ith resampling. 



where N 

The projected 2PCFs u'p(rp) for all the galaxy samples described above are shown in Fig. [201 with different panels 
corresponding to galaxies in different stellar mass bins, as indicated. Open circles with error bars correspond to the 
'Volumel' samples, while solid and dashed lines correspond to the 'Volume2' and 'Mass-limited' samples, respectively. 
Note that the results for the same stellar mass bin are very similar regardless of the sample used, except in the lowest 
stellar mass bin. The big difference between the 'Volume2' from other samples in the lowest stellar mass bin is probably 
produced by the fact that the excess galaxies in other samples are mainly blue galaxies which are relatively bright. 
In general, 'Volume2' sample is unbiased with respect to galaxy color, however suffer more severely from the survey 
volume effect. Therefore, for our model constraints in the main text, we obtain the averaged projected 2PCFs over 
the 'Volumel', 'Volume2' and 'Mass-limited' samples. The resulting covariance matrixes are updated using the total 
600 bootstrap resampling values as well. Thus obtained values of Wp(rp) are listed in Table [5] The related covariance 
matrixes are available from the authors upon request. As we have tested, our results are in good agreement with those 
obtained by Li et al. (2006) from the SDSS DR2 data and those obtained by Guo et al. (2011) from the SDSS DR7 
data, which are based on stellar masses estimated from galaxy spectra and photometries, respectively. For consistency 
with our stellar mass functions, which are all based on the stellar masses obtained using the model of Bell et al. (2003), 
we use our own measurements of u'p(rp) together with the covariance matrixes based on the same stellar masses. 

B. THE STELLAR MASS FUNCTIONS AND CONDITIONAL STELLAR MASS FUNCTIONS 

MEASURED FROM SDSS DR7 

We have constructed group catalougues from the latest SDSS DR7 using the same method described in Yang et al. 
(2007) 0. Here we provide the updated stellar mass functions and conditional stellar mass functions obtained from 
these group catalougues. Readers are referred to Yang et al. (2009b) for the details of these measurements. 

Listed in Table IH] are the stellar mass functions of galaxies. Galaxies are classified into three categories according 
to their memberships in groups: all galaxies (ALL), central galaxies (CENTRAL) and satellites (SATELLITE). For 
galaxies in each category , results are provided for all galaxies (all), red galaxies (red) and blue galaxies (blue), using 
the same color separation as in Yang et al. (2009b). 

For the present paper, investigation in this study, only data for all galaxies in 'ALL' groups are used. This stellar 
mass function is shown in Fig. [2] as open circles with error bars. Fitting the observational data with the Schechter 
function, 

(a+l) 



$(M*) 



exp 



' A4* 



(Bl) 



we obtain the best fit parameters, 0* — 0.0083635, a = —1.117, and log A/* = 10.673. The best fit is shown as the 
solid line in Fig.[5J 

Listed in Table [7] are the conditional stellar mass functions of galaxies. Here halo masses are obtained assuming the 
mass function of the WMAP7 cosmology. 
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TABLE 6 

The galaxy stellar mass functions $(M*) 
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.2612 


0, 


,3709 


± 


0, 


,1511 


1, 


.8474 


± 


0, 


.2058 


1, 


.0873 


it 


0, 


.2604 


0, 


,4294 


it 


0, 


,1365 


0, 


.6578 


it 


0, 


.1831 


8.7 


3, 


.1321 


± 


0, 


.3100 





.8215 


± 


0, 


,1561 


2, 


.3106 


± 


0, 


.2224 


2, 


.1598 


± 


0, 


.2294 


0, 


,3756 


± 


0, 


,1010 


1, 


.7842 


± 


0, 


.1819 


0, 


.9723 


it 


0, 


.1686 


0, 


,4459 


it 


0, 


,1020 


0, 


.5264 


it 


0, 


.1109 


8.8 


3, 


.0391 


± 


0, 


.2499 





.8716 


± 


0, 


,1181 


2, 


.1675 


± 


0, 


.1865 


1, 


.8100 


± 


0, 


.1428 


0, 


,3005 


± 


0, 


,0669 


1, 


.5095 


± 


0, 


.1253 


1, 


.2291 


it 


0, 


.1729 


0, 


,5711 


it 


0, 


,0936 


0, 


.6580 


it 


0, 


.1253 


8.9 


2, 


.7949 


± 


0, 


,2433 





.8404 


± 


0, 


,1442 


1, 


,9545 


± 


0, 


.1538 


1, 


,7266 


± 


0, 


.1265 


0, 


,2997 


± 


0, 


,0582 


1, 


.4269 


± 


0, 


,1027 


1, 


.0683 


it 


0, 


,1745 


0, 


,5407 


it 


0, 


,1141 


0, 


,5276 


it 


0, 


.0929 


9.0 


3, 


,1430 


± 


0, 


,1822 





.9815 


± 


0, 


,1210 


2, 


,1614 


± 


0, 


,1161 


1, 


,9179 


± 


0, 


,0875 


0, 


,3476 


± 


0, 


,0524 


1, 


,5702 


± 


0, 


,0968 


1, 


,2251 


it 


0, 


,1483 


0, 


,6339 


it 


0, 


,0941 


0, 


,5912 


it 


0, 


,0871 


9.1 


3, 


,1047 


± 


0, 


,2357 


1 


.0438 


± 


0, 


,1518 


2, 


,0609 


± 


0, 


,1199 


1, 


,8162 


± 


0, 


,1129 


0, 


,3595 


± 


0, 


,0608 


1, 


,4568 


± 


0, 


,0765 


1, 


,2884 


± 


0, 


,1579 


0, 


,6843 


it 


0, 


,1131 


0, 


,6042 


± 


0, 


,0723 


9.2 


2, 


,9365 


± 


0, 


,1816 


1 


.0557 


± 


0, 


,1398 


1, 


,8808 


± 


0, 


,0760 


1, 


,6895 


± 


0, 


,1021 


0, 


,3411 


± 


0, 


,0622 


1, 


,3484 


± 


0, 


,0606 


1, 


,2470 


± 


0, 


,1151 


0, 


,7146 


it 


0, 


,0969 


0, 


,5324 


it 


0, 


,0409 


9.3 


2, 


,8092 


± 


0, 


,1786 


1 


.0230 


± 


0, 


,1329 


1, 


,7861 


± 


0, 


,0766 


1, 


,5992 


± 


0, 


,0730 


0, 


,3624 


± 


0, 


,0458 


1, 


,2368 


± 


0, 


,0476 


1, 


,2100 


it 


0, 


,1358 


0, 


,6607 


it 


0, 


,1043 


0, 


,5493 


it 


0, 


,0518 


9.4 


2, 


,8013 


± 


0, 


,0925 


1 


.0764 


± 


0, 


,0703 


1, 


,7249 


± 


0, 


,0477 


1, 


,6116 


± 


0, 


,0621 


0, 


,4012 


± 


0, 


,0381 


1, 


,2104 


± 


0, 


,0420 


1, 


,1897 


it 


0, 


,0549 


0, 


,6753 


it 


0, 


,0479 


0, 


,5145 


it 


0, 


,0292 


9.5 


2, 


,5093 


± 


0, 


,1140 


1 


.0816 


± 


0, 


,0917 


1, 


,4277 


± 


0, 


,0418 


1, 


,4360 


± 


0, 


,0522 


0, 


,4290 


± 


0, 


,0420 


1, 


,0070 


± 


0, 


,0354 


1, 


,0733 


it 


0, 


,0804 


0, 


,6526 


it 


0, 


,0637 


0, 


,4208 


it 


0, 


,0290 


9.6 


2, 


,3481 


± 


0, 


,1002 


1 


.0112 




0, 


,0787 


1, 


,3369 


± 


0, 


,0362 


1, 


,3756 


± 


0, 


,0508 


0, 


,4339 


± 


0, 


,0302 


0, 


,9416 


± 


0, 


,0307 


0, 


,9725 


it 


0, 


,0653 


0, 


,5772 


it 


0, 


,0598 


0, 


,3953 


it 


0, 


,0173 


9.7 


2, 


,0970 


± 


0, 


,0640 


1 


.0132 


± 


0, 


,0488 


1, 


,0837 


± 


0, 


,0286 


1, 


,2562 


± 


0, 


,0420 


0, 


,4760 


± 


0, 


,0288 


0, 


,7802 


± 


0, 


,0231 


0, 


,8408 


it 


0, 


,0368 


0, 


,5373 


it 


0, 


,0304 


0, 


,3035 


± 


0, 


,0153 


9.8 


1, 


,9927 


± 


0, 


,0653 


1, 


.0453 




0, 


,0526 


0, 


,9473 


± 


0, 


,0239 


1, 


,2189 


± 


0, 


,0254 


0, 


,5060 


± 


0, 


,0210 


0, 


,7129 


± 


0, 


,0171 


0, 


,7738 


it 


0, 


,0509 


0, 


,5393 


it 


0, 


,0407 


0, 


,2345 


it 


0, 


,0168 


9.9 


1, 


.8551 


± 


0, 


.0555 


1, 


.0426 


± 


0, 


,0446 


0, 


.8125 


± 


0, 


.0210 


1, 


.1423 


± 


0, 


.0240 


0, 


,5284 


± 


0, 


,0176 


0, 


.6139 


± 


0, 


.0143 


0, 


.7128 


it 


0, 


.0398 


0, 


,5142 


it 


0, 


,0337 


0, 


.1986 


it 


0, 


.0119 


10.0 


1, 


.7485 


± 


0, 


.0555 


1, 


.0329 




0, 


,0448 


0, 


.7156 


± 


0, 


.0184 


1, 


.1068 


± 


0, 


.0241 


0, 


,5713 


± 


0, 


,0197 


0, 


.5355 


± 


0, 


.0110 


0, 


.6417 


it 


0, 


.0393 


0, 


,4616 


it 


0, 


,0320 


0, 


.1801 


it 


0, 


.0119 


10.1 


1, 


.6715 


± 


0, 


.0430 


1, 


.0297 


± 


0, 


,0343 


0, 


.6418 


± 


0, 


.0156 


1, 


.0844 


± 


0, 


.0161 


0, 


,5863 


± 


0, 


,0100 


0, 


.4981 


± 


0, 


.0121 


0, 


.5871 


it 


0, 


.0326 


0, 


,4434 


it 


0, 


,0293 


0, 


.1436 


it 


0, 


.0066 


10.2 


1, 


.6340 


± 


0, 


.0417 


1, 


.0560 


± 


0, 


,0331 


0, 


.5780 


± 


0, 


.0142 


1, 


.0963 


± 


0, 


.0122 


0, 


,6465 


± 


0, 


,0088 


0, 


.4498 


± 


0, 


.0082 


0, 


,5377 


it 


0, 


.0346 


0, 


,4095 


it 


0, 


,0289 


0, 


.1282 


it 


0, 


.0087 


10.3 


1, 


,5273 


± 


0, 


,0419 


1 


.0368 


± 


0, 


,0355 


0, 


,4905 


± 


0, 


,0113 


1, 


,0326 


± 


0, 


,0104 


0, 


,6491 


± 


0, 


,0078 


0, 


,3835 


± 


0, 


,0069 


0, 


,4947 


it 


0, 


,0356 


0, 


,3877 


it 


0, 


,0318 


0, 


,1071 


it 


0, 


,0064 


10.4 


1, 


,3308 


± 


0, 


,0339 





.9331 


± 


0, 


,0266 


0, 


,3978 


± 


0, 


,0105 


0, 


,9275 


± 


0, 


,0104 


0, 


,6129 




0, 


,0081 


0, 


,3146 


± 


0, 


,0067 


0, 


,4033 


it 


0, 


,0273 


0, 


,3202 


it 


0, 


,0232 


0, 


,0831 


± 


0, 


,0057 


10.5 


1, 


,0870 


± 


0, 


,0292 





.7817 


± 


0, 


,0237 


0, 


,3052 


± 


0, 


,0084 


0, 


,7882 


± 


0, 


,0095 


0, 


,5383 


± 


0, 


,0066 


0, 


,2499 


± 


0, 


,0051 


0, 


,2988 


it 


0, 


,0225 


0, 


,2435 


it 


0, 


,0193 


0, 


,0553 


it 


0, 


,0045 


10.6 


0, 


,8692 


± 


0, 


,0265 





.6337 


± 


0, 


,0210 


0, 


,2354 


± 


0, 


,0077 


0, 


,6445 


± 


0, 


,0097 


0, 


,4523 


± 


0, 


,0069 


0, 


,1922 


± 


0, 


,0047 


0, 


,2247 


it 


0, 


,0194 


0, 


,1814 


it 


0, 


,0164 


0, 


,0432 


it 


0, 


,0040 


10.7 


0, 


,6629 


± 


0, 


,0208 





.4876 


± 


0, 


,0164 


0, 


,1753 


± 


0, 


,0061 


0, 


,5122 


± 


0, 


,0089 


0, 


,3643 


± 


0, 


,0064 


0, 


,1479 


± 


0, 


,0039 


0, 


,1507 


it 


0, 


,0138 


0, 


,1232 


it 


0, 


,0116 


0, 


,0274 


it 


0, 


,0029 


10.8 


0, 


,4749 


± 


0, 


,0168 





.3555 


± 


0, 


,0133 


0, 


,1194 


± 


0, 


,0045 


0, 


,3796 


± 


0, 


,0078 


0, 


,2776 


± 


0, 


,0056 


0, 


,1020 


± 


0, 


,0031 


0, 


,0953 


it 


0, 


,0101 


0, 


,0779 


it 


0, 


,0088 


0, 


,0173 


it 


0, 


,0019 


10.9 


0, 


,3130 


± 


0, 


,0133 





.2368 


± 


0, 


,0104 


0, 


,0762 


± 


0, 


,0038 


0, 


,2598 


± 


0, 


,0083 


0, 


,1923 


± 


0, 


,0063 


0, 


,0675 




0, 


,0029 


0, 


,0532 


± 


0, 


,0057 


0, 


,0445 


it 


0, 


,0047 


0, 


,0087 


it 


0, 


,0013 


11.0 


0, 


,1913 


± 


0, 


,0086 





.1491 


± 


0, 


,0066 


0, 


,0422 


± 


0, 


,0026 


0, 


,1636 


± 


0, 


,0059 


0, 


,1260 


± 


0, 


,0044 


0, 


,0376 


± 


0, 


,0021 


0, 


,0277 


± 


0, 


,0032 


0, 


,0231 


it 


0, 


,0027 


0, 


,0046 


± 


0, 


,0007 


11.1 


0, 


,1055 


± 


0, 


,0056 


0, 


.0840 


± 


0, 


,0041 


0, 


,0215 


± 


0, 


,0019 


0, 


,0943 


± 


0, 


,0044 


0, 


,0745 


± 


0, 


,0031 


0, 


,0198 


± 


0, 


,0015 


0, 


,0112 


it 


0, 


,0015 


0, 


,0095 


it 


0, 


,0012 


0, 


,0016 


it 


0, 


,0004 


11.2 


0, 


.0540 


± 


0, 


.0028 





.0447 


± 


0, 


,0021 


0, 


.0092 


± 


0, 


.0009 


0, 


.0495 


± 


0, 


.0023 


0, 


,0408 


± 


0, 


,0017 


0, 


.0087 


± 


0, 


.0008 


0, 


.0045 


it 


0, 


.0006 


0, 


,0039 


it 


0, 


,0006 


0, 


.0006 


it 


0, 


.0001 


11.3 


0, 


.0245 


± 


0, 


.0015 





.0207 




0, 


,0013 


0, 


.0039 


± 


0, 


.0004 


0, 


.0231 


± 


0, 


.0014 


0, 


,0194 


± 


0, 


,0011 


0, 


.0037 


± 


0, 


.0004 


0, 


.0015 


it 


0, 


.0002 


0, 


,0013 


it 


0, 


,0002 


0, 


.0001 


it 


0, 


.0001 


11.4 


0, 


.0104 


± 


0, 


.0007 





.0086 




0, 


,0006 


0, 


.0019 


± 


0, 


.0002 


0, 


.0101 


± 


0, 


.0006 


0, 


,0083 


± 


0, 


,0005 


0, 


.0018 


± 


0, 


.0002 


0, 


.0003 


it 


0, 


.0001 


0, 


,0003 


it 


0, 


,0001 


0, 


.0000 


it 


0, 


.0000 


11.5 


0, 


.0042 


± 


0, 


.0003 





.0034 




0, 


,0003 


0, 


.0008 


± 


0, 


.0001 


0, 


.0041 


± 


0, 


.0003 


0, 


,0033 


± 


0, 


,0003 


0, 


.0008 


± 


0, 


.0001 


0, 


.0001 


it 


0, 


.0000 


0, 


,0001 


it 


0, 


,0000 


0, 


.0000 


it 


0, 


.0000 


11.6 


0, 


.0013 


± 


0, 


.0001 





.0010 


± 


0, 


,0001 


0, 


.0003 


± 


0, 


.0001 


0, 


.0013 


± 


0, 


.0001 


0, 


,0010 


± 


0, 


,0001 


0, 


.0003 


± 


0, 


.0001 


0, 


.0000 


it 


0, 


.0000 


0, 


,0000 


it 


0, 


,0000 


0, 


.0000 


it 


0, 


.0000 


11.7 


0, 


.0003 


± 


0, 


,0001 





.0002 


± 


0, 


,0001 


0, 


,0001 


± 


0, 


.0000 


0, 


,0003 


± 


0, 


.0001 


0, 


,0002 


± 


0, 


,0001 


0, 


.0001 


± 


0, 


,0000 


0, 


.0000 


it 


0, 


,0000 


0, 


,0000 


it 


0, 


,0000 


0, 


,0000 


it 


0, 


.0000 



Note. — Column (1): the median of the logarithm of the galaxy stellar mass with bin width A log — 0.05. Column (2 - 4): the stellar mass functions of all, red and blue for 'ALL' group 
members. Column (5-7): the stellar mass functions of all, red and blue for 'CENTRAL' group members. Column (8 - 10): the stellar mass functions of all, red and blue for 'SATELLITE' group 
members. Note that all the galaxy stellar mass functions listed in this table are in units of 10~^/i^Mpc~'^d log M* , where log is the 10 based logarithm. 



TABLE 7 

The conditional stellar mass functions op galaxies 3>(M, IA/^) 



(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


logM* 


[12.1, 12.4] 


[12.4, 12.7] 


[12.7, 13.0] 


[13.0, 13.3] 


[13.3, 13.6] 


[13.6, 13.9] 


[13.9, 14.2] 


[14.2, 14.5] 


> 14.5 



8.9 


0. 


.857 


± 


0, 


,238 


1, 


,534 


± 


0, 


.558 


3 


.944 


± 


1, 


.745 


7, 


.232 


it 


3, 


,599 


0.000 it 0.000 


0.000 it 0.000 


9.0 


0. 


,867 


± 


0, 


,191 


1, 


,793 


± 


0, 


.364 


2 


.569 


± 


0, 


,723 


7, 


,962 


it 


1, 


,834 


10.542 it 4.496 


19.005 it 6.481 


9.1 


0, 


,860 


± 


0, 


,138 


1, 


,117 


± 


0, 


,237 


2 


.984 


± 


0, 


,662 


7, 


,154 


it 


1, 


,573 


12.582 it 4.204 


22.436 it 4.295 


9.2 


0, 


,740 


± 


0, 


,099 


1, 


,353 


± 


0, 


,209 


2 


.221 


it 


0, 


,437 


6, 


,045 


it 


0, 


,845 


11.677 it 3.663 


20.303 ± 4.000 


9.3 


0, 


,626 


± 


0, 


,080 


1, 


,278 


± 


0, 


,180 


2 


.412 


± 


0, 


,321 


4, 


,447 


± 


0, 


,903 


7.831 


± 


1.935 


19.684 ± 3.393 


9.1 


0, 


,572 


± 


0, 


,071 


1, 


,252 


± 


0, 


,136 


3 


.118 




0, 


,316 


4, 


,153 


± 


0, 


,540 


8.288 


± 


1.085 


20.259 ± 2.992 


9.5 


0, 


,476 


± 


0, 


,057 


1, 


,064 


± 


0, 


,110 


2 


.034 


± 


0, 


,222 


4, 


,965 


± 


0, 


,653 


9.177 




1.263 


21.391 it 2.823 


9.6 


0, 


,522 


± 


0, 


,054 


0, 


,991 


± 


0, 


,093 


2 


.512 


± 


0, 


,218 


4, 


,064 


± 


0, 


,455 


8.755 




1.012 


16.330 it 1.757 


9.7 


0, 


,414 


± 


0, 


,037 


0, 


,873 




0, 


,078 


1 


.691 


± 


0, 


,160 


3, 


,553 


± 


0, 


,382 


6.781 




1.065 


15.900 ± 1.777 


9.8 


0, 


,419 




0, 


,046 


0, 


,817 


± 


0, 


,088 


1 


.700 


± 


0, 


,1.38 


3, 


,288 


± 


0, 


,283 


7.694 




0.865 


14.748 ± 1.399 


9.9 


0, 


,401 


± 


0, 


,028 


0, 


,802 


± 


0, 


,055 


1, 


.733 




0, 


,133 


2, 


,959 


it 


0, 


,243 


7.274 


± 


1.067 


16.538 it 1.313 


10.0 


0. 


,339 


± 


0, 


,052 


0, 


,766 


± 


0, 


.057 


1, 


.598 


± 


0, 


.104 


2, 


.983 


it 


0, 


,197 


6.379 


± 


0.630 


15.380 it 1.293 


10.1 


0. 


,283 


± 


0, 


,055 


0, 


,673 


± 


0, 


.037 


1, 


.388 


± 


0, 


.073 


3, 


.167 


it 


0, 


,173 


6.209 


± 


0.453 


13.795 it 0.807 


10.2 


0. 


,176 


± 


0, 


,040 


0, 


,610 


± 


0, 


.039 


1, 


.391 




0, 


.079 


2, 


.850 


it 


0, 


,133 


5.491 


± 


0.486 


13.024 it 0.811 


10.3 


0, 


,088 


± 


0, 


,034 


0, 


,547 


± 


0, 


.098 


1, 


.300 


± 


0, 


.097 


2, 


.473 


it 


0, 


,147 


5.500 


± 


0.242 


11.631 it 0.710 


10.4 


0, 


,010 


± 


0, 


,013 


0, 


,357 


± 


0, 


,061 


1 


.077 


± 


0, 


,071 


2, 


,332 


it 


0, 


,130 


4.577 




0.218 


10.016 it 0.700 


10.5 


0, 


,002 


± 


0, 


,002 


0, 


,113 


± 


0, 


,039 





.828 


± 


0, 


,120 


1, 


,918 
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Note. — Column (1): the median of the logarithm of the galaxy stellar mass with bin width A log M* — 0.05. Column (2 - 10): the conditional stellar mass functions in 
halos of different mass ranges as indicated. The average halo masses in these bins are log < >— 12.26, 12.55, 12.85, 13.15, 13.44, 13.74, 14.04, 14.33, 14.72, respectively. Note 
that all the conditioanl galaxy stellar mass functions listed in this table are in units of dlogM*, where log is the 10 based logarithm. Here results are listed for satellite (upper 
part) and central (lower part) galaxies separately, which can be distinguished with zero measurements. Note here the stellar masses for central galaxies should be converted using 
logM,,c = logM, - 1.0. 
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